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ABSTRACT 

We report on variations in important X-ray emission lines in a series of Chandra grating spectra 
of the supermassive colliding wind binary star rj Car, including key phases around the X-ray min- 
imum/periastron passage in 2003.5. The X-rays arise from the collision of the slow, dense wind of 
rj Car with the fast, low-density wind of an otherwise hidden companion star. The X-ray emission 
lines provide the only direct measure of the flow dynamics of the companion's wind along the wind- 
wind collision zone. We concentrate here on the silicon and sulfur lines, which are the strongest and 
best resolved lines in the X-ray spectra. Most of the line profiles can be adequately fit with sym- 
metric Gaussians with little significant skewness. Both the silicon and sulfur lines show significant 
velocity shifts and correlated increases in line widths through the observations. The TZ — forbidden- 
to-intercombination ratio from the Si XIII and S XV triplets is near or above the low-density limit 
in all observations, suggesting that the line- forming region is > 1.6 stellar radii from the companion 
star. We show that simple geometrical models cannot simultaneously fit both the observed centroid 
variations and changes in line width as a function of phase. We show that the observed profiles can be 
fitted with synthetic profiles with a reasonable model of the emissivity along the wind-wind collision 
boundary. We use this analysis to help constrain the line formation region as a function of orbital 
phase, and the orbital geometry. 

Subject headings: X-rays: stars -stars: early-type-stars: individual {rj Car) 



1. INTRODUCTION 



Th e supermassive star 77 Car ([Davidson fc Humphreys! 
Il997f ) is notorious for its extraordinarily large luminos- 
ity and its im plicitly large mass (L > 4 x IQ^Lq and 
M - lOOMg. iHiriier et al.l [200lh . the beautif ul bipo- 
lar " Homunculus" nebula which shrouds it (|Gaviolal 
|195C1| ). its wild instability (most notably the "Great 
Eruption" of 1843 which created th e Homunculus) and 
its continued broad-b and variations ijSterken et al.lll996l : 
[Davidson et al.l [19991 ). Understanding 77 Car is impor- 
tant for a wide variety of astrophysical topics regarding 
the formation and evolution of extremely massive stars, 
the processes by which such stars lose mass and angular 
momentum, and the ways in which they interact with 
their surroundings. 

rj Car exhibits variability over a wide rang e of wave- 
lengths, from radio (Duncan fc White! l2003l). th rough 
infrared (IWhitelock et all 11994 l200l fP aminch Il996t 
iDamineh et al." 1997*. '20001: iDavidson et al.il200a), optical 
JStcincr & Damincli 2 0(3^, and u ltravi olet (ISmith et aTl 
[2004) to X-rays (Ishibashi et al.l 119991 : ICorcoranll20'05i r 
All these variations have a characteristic cycle of almost 
exactly 2024 days, which strongly sugg ests that r? Car 
is a long period (P = 2024 day) binary (|Daminelilll996t 
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iDamineli et al.! |1997[ ) . The observed variability is be- 
lieved to result (directly or indirectly) from the interac- 
tion of the fast wind (fc 3000 km s~ ) and ionizing 
radiation from the companion with the dense, slow wind 
of the Luminous Blue Variable (LBV) primary (wj, ~ 
500 km s~^). In this scheme, the X-rays are produced 
by the collision of the two stars' winds, which causes 
the c ompanion's fast wind to be shock-heated to t ens of 
MK (iPittard et aLlll998HPittard fc Corcor"Mll2002[ ). The 
high temperature of the shocked wind of the companion 
explains the hard X-rays {kT > 4 k eV) first directly as - 
sociated with the star by Einstein (jSeward "eral][l979h . 
Similar hard X-ray emission is seen from WR 140, the 
"canoni cal" long period ecc entric massive colliding wind 
binary (jPollock et al.ll2005[ ). 

Our understanding of the system has become more so- 
phisticated due in part to dense multiwavelength mon- 
itoring near the times of the X-ray minima in 1998 
and 2003.5. These observations showed that, at the 
same time that the X-ray brightness of the source 
reaches minimum, the ionization state of the circum- 
stellar medium rapidly decrease s (iDuncan et al.l |199 
iNielsen et aLll2007[ ). t he infrared (IWhitelock et al. 200'_, 
and millimeter- wave ([Abraham et al.ll2005bl ) brightness 
of the source also drops, absorption components in ex- 
cited He I P-Cygni e mission lines undergo rapid blue- 
to-red vel ocitv shifts (iNielsen et al.l[2007i)" He II 4686- A 
emission (jSteiner fc Damineli! 12001 iMartin et al.1 12006!) 
appears, shows a similar blue-to-red centroid shift, then 
disappears, and th e far UV flux from rj Car drops rapidly 
(jlping et all 120051 ). In X -rays, the hottest electron tem- 
perature stayed the same, b ut the ionization balanc e of 
Fe ions changed remarkably (jHamaguchi et aLll2007f l. In 
all colliding wind models these changes (which last only 
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about 90 days of the 2024-day cycle) occur near perias- 
tron passage, and require a high eccentricity (e ~ 0.9). 
However, important details regarding the nature of the 
wind-wind collision are still not well constrained; there 
is still debate concerning, for example, whether the X- 
ray minimum occurs near inferior conjunction (when the 
companion is in front of the LBV primary) or superior 
conjunction; the magnitude of the companion's wind ve- 
locity; and the mass loss rates from either star. These 
uncertainties limit our understanding of how the com- 
panion star affects the system, and, ultimately limit our 
knowledge of the evolutionary state of the system. 

The detailed analysis of excited He I P- Cygni absorp- 
tion li nes in spatially resolved spectra by iNielsen et al.l 
I2OO7') showed radial velocity variations which mimic the 
orbital radial velocity variations expected in an eccentric 
(e « 0.9) binary system with the semi-major axis pointed 
towards the observer (longitude of periastron uj ~ 270°) 
and an assumed inclination i k, 41°. These spectral vari- 
ations suggest that the ionized helium zone in the wind 
of the cool, massive primary star approaches the observer 
prior to periastron passage. They also showed that the 
velocity amplitude of the He I P-Cygni absorption com- 
ponents (^^140 km s^^) was larger than expected if the 
absorption arises in the dense wind of the more massive 
star. They concluded that the velocity variations are 
probably strongly influenced by ionization effects due to 
the interaction of the companion star's photospheric UV 
radiation with the wind of the cool primary star. They 
also suggested that some of the He I emission might orig- 
inate within or near the wind-wind collision and thus 
could be a diagnostic of that collision. However the com- 
plex influence of the companion's radiation with the pri- 
mary wind makes interpretation of such diagnostics far 
from straightforward. (For an alternative explanation of 
the He I observations, in which the lines are assumed 
to form in the acce leration zone of the secondary, see 
iKashi k SokiJ[2007l l 

X-ray line profiles provide the most direct probe of 
the dynamics of the wind of the unseen companion af- 
ter it is shock-heated in the wind-wind interaction, since 
these lines originate in the high temperature plasma near 
the wind-wind shock interface. X-ray lines directly re- 
flect the dynamic properties of this hot shocked gas. In 
this paper we present our analysis of the high resolu- 
tion X-ray grating spectra of 77 Car obtained by the High 
Energy Transmissio n Grating Spectrometer (HETGS; 
iMarkert et al.lll994D on the Chandra X-ray Observatory 
( Weisskopf et al.l |2002| ) obtained as part of a large ob- 
serving campaign around the time of the 2003.5 X-ray 
minimum . A pre li mina ry analysis of these data has ap- 
peared in lHenlevI (|2005f ). 

In this paper we discuss our analysis of spectra in 
the energy range near 2 keV obtained by the Medium 
and High Energy Gratings (MEG and HEG). This en- 
ergy range is dominated by line emission from Si and 
S hydrogen- like and helium- like ions. These lines form 
in the cooler regions of the shocked gas farther along 
the wind-wind collision zone, and thus provide a better 
measure of the flow dynamics of the shock-heated wind 
of the companion along the colliding wind interface than 
the iron lines, which originate near the hottest part of the 
shock near the stagnation point where flow velocities arc 



low. In this energy range the HETGS first order spectra 
has sufficient resolution to resolve the component lines of 
the He-like triplets providing useful density and temper- 
ature diagnostics. Unfortunately, potentially crucial line 
emission from C, N, and O (which could be used to mea- 
sure abundances of the shocked companion's wind and 
help constrain the evolutionary state of the companion) 
are not observable in the central source due to the heavy 
absorption by the cold gas and dust in the Homunculus. 

This paper is organized as follows. The observations 
and the data reduction are described in fJS] and the 
HETGS silicon and sulfur emission lines are discussed 
in SjHl In f|4] we apply a simple geometrical model of the 
wind-wind collision to the variations in line centroids and 
widths. In Sj5]we apply synthetic colliding wind line pro- 
files to the observed HETGS silicon and sulfur profiles. 
We discuss the results of the emission line analysis in SJHI 
and our conclusions are presented in [J7l Throughout this 
paper we quote la errors. 

2. OBSERVATION DETAILS AND DATA REDUCTION 

The details of the six Chandra HETGS observations of 
77 Car are given in Table [T] For the purposes of this pa- 
per we designate the observations with CXO, subscripted 
with the date in YYMMDD format (cf. iHamaguchi et aH 
l2007t l. The earliest observation was in 2000 November 
(CXO noiiig; ICorcoran et"al]|2001bl: iPittard fc CorcoranI 
l2002l ). approximately half-way between the previous X- 
ray minimum in late 1997 and the X-ray minimum in 
mid-2003. The second observation was taken approx- 
imately one year later (2002 October; CXOo2ioi6)j by 
which time the X-ray flux had increased by a factor 
of '^2. The four remaining observations were taken 
over the space of approximately five months around 
the X-ray minimum which occurred, as expected, in 
late 2003 June. In particular, they approximately cor- 
respond to X-ray maximum (2003 May; CXO030502), 
the early part of the descent to X-ray minimum (2003 
June; CXOo306i6)j the X-ray minimum itself (2003 July; 
CXOo3072o)i a-nd the recovery from the minimum (2003 
September; CXOo30926)- AH data were read out using 
the Advanced Camera for Imaging Spectroscopy spec- 
troscopic array (ACIS-S). The outer ACIS-S CCD chips 
(SO and S5) were switched off, and we used a reduced 
read-out window in order to reduce pileup. This trun- 
cates the low-energy spectra but results in little real data 
loss since the stellar source is heavily absorbed. The 
spectra at energies E <?> keV obtained during and just 
after the X-ray minimum (CXO030720 and CXO030926) 
are contaminated by the "Central Const ant Eni i ssion" 
(CCE) component identified bv fHamaguchi et a l.l (j2007D 
from XMM-Newton observations taken during the 2003 
X-ray minimum. This means that the silicon and sulfur 
lines from these two spectra do not accurately reflect the 
emission from the colliding wind plasma alone (with the 
exception of S xvi in CXO030926: which is not as badly 
contaminated). However, for completeness, we include 
measurements of the line properties for all six observa- 
tions, including CXO030720 and CXO0309261 in our dis- 
cussion in §||3]and[4l 

The data for all six observations were reduced from 
the Level 1 events files using CIAO"^ v3.4 and CALDB 

http://cxc.harvard.edu/ciao 
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TABLE 1 

Chandra-RET GS Observations of -q Carinas 



Observation 


Observation 


Start 


Phase'^ 


Exposure 




HEG 




MEG 




IDb 


date 




(ks) 


Counts^ 


Rate (s^i) 


Counts'^ 


Rate (s-i) 


CXO001119 


632 


2000 Nov 19 


0.528 


89.5 


18459 


0.206 


20772 


0.232 


CXO021016 


3749 


2002 Oct 16 


0.872 


91.2 


38160 


0.418 


45038 


0.493 


CXOo30502 


3745 


2003 May 2 


0.970 


94.5 


78264 


0.828 


81925 


0.867 


CXO0306I6 


3748 


2003 Jun 16 


0.992 


97.2 


42411 


0.436 


40553 


0.417 


CXOo30720 


3746 


2003 Jul 20 


1.009 


90.3 


1183 


0.013 


1725 


0.019 


CXOo30926 


3747 


2003 Sep 26 


1.043 


70.1 


11137 


0.159 


8098 


0.116 



^ Observation identification used in this paper Rafter [Hamaguchi et 3111200711 . 



Official Chandra observation identification. 
^ Mid-observation phase, calculated using the emphemeris in lCorcora 
^ Total number of first-order ( + 1 and —1) non-background-subtraeted eounts. 



vS.S.O.l. These versions are much impr o ved ove r the 
earlier versions used by ' Corcoran et al.l ()2001bf ) and 
iHenlevI {2005). We followed the threads available 
from the Chandra website^. We first removed the 
acis_d.etect_af terglow correction, and generated a 
new bad pixel file using acisjrunjiotpix. We then 
reprocessed the Level 1 events file with the latest cal- 
ibration using acis_process_events. This applies a 
new ACIS gain map, the time-dependent ACIS gain 
correction, the ACIS charge transfer inefficiency (CTI) 
correction, and pixel and PHA randomization. We 
then used tgdetect to determine the position of the 
zeroth-ordcr image of 77 Car, tg_create_mask to deter- 
mine the location of the HEG and MEG "arms" , and 
tgjresolve_events to assign the measured events to the 
different spectral orders. After applying grade filters 
{ASCA grades 0, 2, 3, 4, and 6 were kept) and good time 
intervals, we used destreak to remove streaks caused by 
a flaw in the serial readout which randomly deposits sig- 
nificant amounts of charge along the pixel row as charge 
is read out. Finally, we used tgextract to extract the 
grating spectra from the events file. Spectral response 
files were also generated following the Chandra threads: 
we generated redistribution matrix files (RMFs) and an- 
cillary response files (ARFs) using mkgrmf and f ullgarf , 
respectively. 

The Chandra HETGS spectra of 77 Car for each of the 
six Chandra observations are shown in Figures [1] through 
[51 For each observation, the +1 and —1 orders of each 
grating (HEG and MEG) have been co-added, and the 
spectra have been binned up to 0.01 A, except for the ob- 
servation taking during the X-ray minimum (CXO030720; 
Fig. [5]), which has been binned up to 0.02 A. Note that 
the spectra are shown with the same y-axis range, except 

for CXOo30720- 

With the exception of CXO030720, which is the faintest 
spectrum by an order of magnitude, the spectra all 
exhibit prominent continuum emission and numerous 
emission lines. Particularly prominent are forbidden- 
intercombination-resonance (f-i-r) triplets from He-like 
Fe XXV (Aresonance = 1-85 A), S XV (Ajos = 5.04 A) and 
Si XIII (Arcs — 6.65 A), Lya emission from H-like S xvi 
(A = 4.73 A) and Si xiv (A = 6.18 A), and K-shell fluo- 
rescent emission from cool Fe (A = 1.94 A). Other lines 
which are visible (not necessarily in all spectra) include 

® http:/ /cxc. harvard.edu/ciao/threads 
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Fig. 1. — Chandra HETGS spectra of rj Car from CXO001119. 
For each grating (HEG and MEG) the -|-1 and —1 orders have been 
co-added, and the spectra have been binned up to 0.01 A. 
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Fig. 2.— As Figure[T] but for CXOo2ioi6. 



4 



HENLEY ET AL. 



CX0„3„,„g 0= 0.970 





0.5 




0.4 




3 


Q 


0.2 














0.1 


c 




< 





o 




CD 






5 








0.4 






o 


0.3 




0.2 




0.1 










° Ca XX-| 


pCa Xixl ' ' ' 


1 1 1 1 1 1 1 _ 








HEG -1/+1 J 
: 


u 






III : 

'Jk'wj^,... .. 1 : 




2 


4 6 


8 10 




2 4 6 8 

Wave-length (Arifrstrom) 

Fig. 3.— As Figure [H but for CXOo30502. 
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Fig. 4. — As Figure [1] but for CXOosoeie- 

Ca XX at 3.0 A, Ca xix at 3.2 A, Ar xviii at 3.7 A, 
Ar XVII + S XVI at 4.0 A, Si xiv at 5.2 A, and Si xiii 
at 5.7 A. The Fe K lines will be discussed elsewhere 
(Paper II; M. F. Corcoran et al., in preparation). Here 
we concentrate on the brightest of the lower-excitation 
lines: the H-like Lya line and He-like f-i-r triplet from Si 
and S. Although line shifts and widths can be measured 
for some of the other lines, the four lines that we con- 
centrate on here are the only ones for which results can 
be obtained from all six observations. Furthermore, the 
analysis of these weaker lines is consistent with the analy- 
sis of the stronger lines presented h ere (for more detailed 
discussion of these weaker lines see lHenlevl[200l . 

3. SILICON AND SULFUR LINE PROFILES 
3.1. Gaussian Modeling 
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Fig. 5. — As Figure [T] but for CXOo30720i with the spectrum 
binned up to 0.02 A, instead of 0.01 A. Note that the ly-axis range 
is different from Figures [T] through |4] and |6] 
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Fig. 6.— As Figure [H but for CXOo30926- 

We analyzed the Chandra spectra of r] Car using un- 
binned, non-co-added spectra, so no spectral informa- 
tion was lost. Because som e bins conta in low numbers of 
counts, the Cash statistic (Cash"1979) was used instead 
of the statistic. To measure each emission line's prop- 
erties, we analyzed each line (or multiplet) individually 
over a narrow range of wavelengths encompassing just 
the line of interest. We then fit a model to the data com- 
prising a power-law continuum component plus Gaussian 
components to model the line emission. The number of 
Gaussians used, and how their parameters are tied to- 
gether, dB2ended_on the nata ana- 
lyzed (jPollock et all 120051 iHmlev et al.ll2005[ ). For the 
Lya lines (which are closely spaced doublets, separated 
by w5 mA), we used two Gaussians. The Doppler shifts 
of the two components were constrained to be equal. 
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as were their widths, and the intensity of the longer- 
wavelength component was constrained to be half that 
of the shorter- wavelength component. The He- like f-i-r 
triplets were fit with three Gaussians, the Doppler shifts 
and widths of which were tied together as for the Lyman 
lines, but the amplitudes of which were allowed to vary. 
For the intercombination line we used the rest wavelength 
of the 2'^Pi — > 1 ^So transition, and ignored the fainter 
2 3P2 ^ 1 iSo transition. 

The analysis described here was carried out using 
SHERPA, as distributed with CIAO v3.4. The data were 
not background subtracted, as the Cash statistic can- 
not be used on background-subtracted data, nor was the 
background separately modeled out. This is not a prob- 
lem because for the lines of interest the background count 
rate is more than an order of magnitude lower than the 
source count rate in the relevant energy range. Further- 
more, the background spectra show no prominent spec- 
tral features, so any background contribution would be 
included in the continuum component used in the fitting. 

Our procedure for a given line from a given observation 
was to fit the same model to all four spectra (HEG ±1, 
MEG ±1) simultaneously. We then assessed goodness-of- 
fit using a Monte Carlo method (as the Cash statistic by 
itself gives no goodness of fit information) , using a similar 
method to that of .Helsdon fc Ponman ( 2000). The best- 
fit model was used to simulate 1000 synthetic emission 
lines. Poisson noise was added to each simulated line, 
and then each was compared with the original model to 
calculate its Cash statistic. Hence, for a given emission 
line model, we obtained the distribution of Cash statistic 
values expected for datasets generated from that model. 
By comparing the observed Cash statistic with this dis- 
tribution, we determined the probability that the model 
could have produced the observed data. In practice we 
did this by measuring the mean and standard deviation 
(cr) of the simulated Cash statistic values - if the ob- 
served Cash statistic lay more than 2a away from the 
mean, the fit was regarded as "poor" . 

We found that, when fitting to all four spectra simulta- 
neously, Gaussian profiles gave acceptable fits to most of 
the lines. A visual inspection of the poorer fits indicated 
that the lines in different orders were sometimes slightly 
offset from each other in wavelength. This may be due to 
uncertainty in the determination of the centroid position 
of the zeroth-order image on the ACIS-S detector - if the 
determined position were offset from the true position, 
the wavelengths in the -1-1 and —1 orders would be offset 
in opposite directions. To overcome this, where possi- 
ble we fit the model to the four spectra individually, and 
then averaged the results. For some fainter lines (the sul- 
fur lines in CXOooiiig, and the lines in CXO030720 and 
CXOo30926) we were unable to constrain the model in 
all four individual spectra. In these cases, we adopted 
the results obtained by fitting all four spectra simultane- 
ously. For the Si xiii triplet in CXO030720, even this did 
not work, and instead we obtained our results by fitting 
just to the MEG -1-1 and —1 spectra. 

The emission line shifts, widths, fluxes, and equiva- 
lent widths measured from this Gaussian modeling are 
given in Table [2l The rest wavelengths are adopted 
from ATOMDB^ vl.3.1. Table [3] shows the results in 

^ http://cxc.harvard.edu/atomdb 



Table [2] expressed as velocities. Figures [7] and |8] show the 
Si XIII and Si xiv lines from the four brightest spectra 

(CXO001119, CXO0210161 CXOo30502, and CXO0306I6), 
along with the best-fitting Gaussian line model. The 
models were fit to each spectral order individually, which 
is why in several cases the Gaussians are offset in the dif- 
ferent orders. 

Figure [9] shows the measured line shifts and widths 
plotted against phase (see Table [ij, where (j) = 1 cor- 
respo nds to the start of the X-ray minimum in 2003 
June (jCorcoranI |2005( ) . We have not corrected for the 
systemic velocity of rj Car (—8 km s~^; ISmithI [20041 ) . 
as it is negligible compared with the measurement er- 
rors. The general trend of the line shifts is that 
the lines have small blueshifts of ^100 km s~^ away 
from the X-ray minimum (CXO001119 and CXO021016)) 
the blueshifts increase to ~300-700 km s~^ just be- 
fore the X-ray minimum (CXO030502 and CXO0306I6; 
note that the lines in CXO0306I6 are generally more 
blueshifted than in CXO030502), and the blueshifts re- 
turn to '^lOO km s^ after the start of X-ray minimum 
(CXOo30720 and CXOo30926)- The exception to this is 
S XVI which is slightly redshifted in the last two obser- 
vations. The general trend of the line widths is that they 
mcrease from ~800 km s"^ (FWHM) away from the X- 
ray minimum to ~1400 km s~ just before the start of 
minimum, and then return to ^800 km s~^ afterward. 
However, as noted in fj2l the last two observations (be- 
ing much fainter than the previous ones) are contami- 
nated by the CCE component (Hamaguchi et al. 200"^ 
at wavelengths longward of about 4A. As a result of this 
contamination, the shifts and widths determined from 
the silicon and sulfur lines in these two spectra do not 
accurately reflect the kinematics of the wind-wind colli- 
sion (with the exception of S xvi in CXOo30926j which is 
not as badly contaminated). 

The connection between the variation in the line shifts 
and the variation in the line widths is further illustrated 
in Figure [101 There is a clear correlation between shift 
and width, with the broader lines being more blueshifted. 
For these data. Spearman's rank correlatio n coefficient 
is —0 .58, and Kendall's r statistic is —0.47 ([Press et al.l 
II992I ). Both of these statistics show that correlation is 
significant at the 1% level. 

Figure [Tl] shows the variation in the emission line 
fluxes, plotted with the 2-10 keV RXTE hghtcurve 
(|Corcoranll2005f l. For S XV and Si Xiii we plot the res- 
onance line flux. As expected, the variation in the line 
fluxes generally follows that of the broadband emission. 
However, not all the lines' fluxes vary in the same way - 
for example, the Si xiii flux does not rise as much as the 
Si XIV flux in the first three observations, which in turn 
does not rise as much as the S xvi flux. These differences 
between the lines are shown more clearly in Figure I12|, 
which shows the ratios of the line fluxes to the contempo- 
raneous 2-10 keV flux measured with RXTE (|C orcorani 
120051 ). The ratios are normalized to the values from 
CXOooiii9- From CXO001119 to CXO021016 (0 = 0.528 
to 0.872), the emission lines stay fairly constant with re- 
spect to the broadband flux (the S xvi and Si xiv lines 
actually brighten slightly). However, just before the X- 
ray minimum (CXO030502 and CXO0306I6; = 0.970 and 
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TABLE 2 
Emission Line Fit Results 



Ion 


Line 


Ao 


Aobs 


AA (FWHM) 


Flux 


EW 


Fitting 






(A) 


(A) 


(mA) 


(10-5 ph cm-2 s-l) 


(A) 


metiiod 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


CXOooiii!) 


S XVI 


Lya 


4.7274 


4.7267 ±0.0012 


20.2 ±4.9 


2.96 ±0.34 


0.041 ±0.005 


(a) 


S XV 


r 


5.0387 


5.0366 ± 0.0007 


13.7±1.8 


4.99 ±0.41 


0.087 ±0.007 


(a) 




i 


5.0665 


5.0644 


13.8 


0.97 ±0.27 


0.017 ±0.005 


(a) 




f 


5.1015 


5.0994 


13.9 


2.41 ±0.33 


0.042 ± 0.006 


(a) 


Si XIV 


Lya 


6.1804 


6.1797 ±0.0006 


16.5 ±1.7 


2.43 ±0.13 


0.117 ±0.006 


(b) 


Si XIII 


r 


6.6479 


6.6461 ± 0.0005 


12.3±1.3 


2.57±0.16 


0.171 ±0.011 


(b) 




i 


6.6882 


6.6864 


12.4 


0.35 ±0.09 


0.024 ± 0.006 


(h) 




r 


().7i():-i 


(i.7:-!8.j 


12. .J 


1.7."') ±0.13 


0.121 ±0.00!) 


(b) 


CXO021016 


S XVI 


Lya 


4.7274 


4.7262 ±0.0006 


14.1 ±2.0 


6.97 ±0.44 


0.042 ± 0.003 


(b) 


S XV 


r 


5.0387 


5.0374 ± 0.0006 


13.6 ±1.5 


9.27±0.67 


0.068 ±0.005 


(b) 




i 


5.0665 


5.0652 


13.7 


3.08 ±0.49 


0.022 ±0.004 


(b) 




f 


5.1015 


5.1002 


13.8 


3.79 ±0.53 


0.028 ±0.004 


(b) 


Si XIV 


Lya 


6.1804 


6.1772 ±0.0005 


23.3 ± 1.6 


5.62 ±0.23 


0.103 ±0.004 


(b) 


Si XIII 


r 


6.6479 


6.6458 ±0.0004 


13.8 ± 1.1 


4.89 ±0.23 


0.133 ±0.006 


(b) 




i 


6.6882 


6.6861 


13.9 


0.50 ±0.12 


0.014 ±0.003 


(b) 




f 


6.7403 


6.7382 


14.0 


3.08 ±0.18 


0.087 ±0.005 


(b) 


CXO()3()5()2 


S XVI 


Lyo 


4.7274 


4.7221 ±0.0006 


21.0 ±1.9 


11.38 ±0.58 


0.043 ±0.002 


(b) 


S XV 


r 


5.0387 


5.0335 ± 0.0009 


22.1 ±2.0 


11.45 ±0.78 


0.066 ±0.005 


(b) 




i 


5.0665 


5.0613 


22.2 


3.34 ±0.59 


0.017 ±0.004 


(b) 




f 


5.1015 


5.0962 


22.4 


5.71 ±0.68 


0.027 ±0.004 


(b) 


Si XIV 


Lya 


6.1804 


6.1698 ±0.0005 


33.7 ±1.4 


7.68 ±0.26 


0.125 ±0.004 


(b) 


Si XIII 


r 


6.6479 


6.6397 ± 0.0006 


24.6 ±1.6 


5.55 ±0.27 


0.122 ±0.005 


(b) 




i 


6.6882 


6.6800 


24.7 


0.62 ±0.16 


0.014 ±0.004 


(b) 




f 


6.7403 


6.7322 


24.9 


3.52 ±0.22 


0.080 ± 0.005 


(b) 


CXO0306I6 


S XVI 


Lya 


4.7274 


4.7213 ±0.0009 


22.6 ±2.8 


4.96 ±0.39 


0.051 ±0.004 


(b) 


S XV 


r 


5.0387 


5.0312 ±0.0008 


21.5±2.2 


8.18 ±0.56 


0.113 ±0.008 


(b) 




i 


5.0665 


5.0590 


21.6 


2.50 ±0.41 


0.036 ± 0.006 


(b) 




f 


5.1015 


5.0939 


21.8 


4.63 ± 0.45 


0.068 ±0.007 


(b) 


Si XIV 


Lya 


6.1804 


6.1667 ±0.0007 


28.3 ± 1.9 


2.92 ±0.15 


0.130 ±0.007 


(b) 


Si XIII 


r 


6.6479 


6.6329 ±0.0009 


35.7 ±2.2 


3.42 ±0.20 


0.244 ±0.014 


(b) 




i 


6.6882 


6.6731 


35.9 


0.34 ±0.13 


0.025 ±0.010 


(b) 




f 


6.7403 


6.7251 


36.2 


2.30 ±0.16 


0.172 ±0.012 


(b) 


CXOo30720 


S XVI 


Lya 


4.7274 


4.7299 ±0.0041 


21.8 ±9.7 


0.19 ±0.06 


0.085 ±0.027 


(a) 


S XV 


r 


5.0387 


5.0374 ±0.0015 


13.8 ±4.0 


0.63 ±0.13 


0.177±0.037 


(a) 




i 


5.0665 


5.0652 


13.9 


0.14 ±0.08 


0.038 ±0.022 


(a) 




f 


5.1015 


5.1002 


14.0 


0.45 ±0.11 


0.113±0.029 


(a) 


Si XIV 


Lya 


6.1804 


6.1766 ±0.0019 


12.1 ±5.3 


0.16 ±0.03 


0.098 ±0.020 


(a) 


Si XIII 


r 


6.6479 


6.6421 ±0.0013 


19.8 ±3.6 


0.58 ±0.07 


0.372 ±0.048 


(c) 




i 


6.6882 


6.6824 


19.9 


0.05 ±0.03 


0.035 ±0.022 


(c) 




f 


6.7403 


6.7344 


20.1 


0.36 ±0.06 


0.248 ±0.039 


(c) 


CXOo30926 


S XVI 


Lya 


4.7274 


4.7291 ±0.0020 


8.0 ±7.0 


0.50±0.11 


0.048 ±0.011 


(a) 


S XV 


r 


5.0.387 


5.0370 ± 0.0012 


11.9 ±3.4 


1.09±0.19 


0.167 ±0.030 


(a) 




i 


5.0665 


5.0648 


12.0 


0.48 ±0.15 


0.076 ±0.023 


(a) 




f 


5.1015 


5.0998 


12.0 


0.48 ±0.14 


0.077 ±0.023 


(a) 


Si XIV 


Lya 


6.1804 


6.1763 ±0.0032 


25 ±10 


0.24 ±0.06 


0.103 ±0.025 


(a) 


Si XIII 


r 


6.6479 


6.6452 ±0.0021 


20.2 ±4.6 


0.25 ±0.06 


0.147 ±0.035 


(a) 




i 


6.6882 


6.6855 


20.3 


0.03 ±0.03 


0.016 ±0.019 


(a) 




f 


6.7403 


6.7376 


20.5 


0.21 ±0.05 


0.128 ±0.032 


(a) 



Note. — Values without quoted errors were tied to other fit parameters. Col. (3): Rest wavelengths from 
ATOMDB vl.3.1. For each Lya line we give the wavelengths of the brighter component; the wavelengths of the 
fainter components are 4.7328 A (S xvi) and 6.1858 A (Si xiv). Col. (4): Observed wavelength. Col. (5): Observed 
line width. Col. (6): Observed line flux. For each Lya line we give the flux of the brighter component; the fluxes of 
the fainter components are half of these values. Col. (7); Equivalent width. Col. (8): (a) Fitting to HEC ±1 and 
MEG ±1 simultaneously; (b) Fitting to HEG ±1 and MEG ±1 individually and averaging the results, (c) Fitting 
to MEG ±1 simultaneously. 
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CX0„3„,„3 Si XIII f-i-r 



cxo„ 



Si XIIT f-i-r 



HEG -1 - 




6.6 6.7 6.8 

Wavelength (Angstrom) 



HEG +1 - 




6.6 6.7 6.B 

Wavelength (Angstrom) 



Fig. 7. — Si xiii f-i-r triplets from CXO001119, CXO0210161 CXO030B02, and CXO0306I61 showing the best-fitting Gaussian line model ob- 
tained by fitting to each spectral order individually. The vertical dashed lines show the rest wavelengths of the resonance, intercombination, 
and forbidden lines (6.6479, 6.6882, and 6.7403 A, respectively). 
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cxo„ 



Si XIV Lya 
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o.ae 
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0.05 
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Si XIV Lya 



CXO„ 



Si XIV Lya 
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0.04 



„ 
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" 
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O 
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0,08 
0.06 
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0.02 








1 1 1 1 1 1 1 1 1 _ 
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faitiiiijiiii^iMi^iijiiii' 




MEG -1 - 




MEG +1 -\ 



6.15 6.2 6,25 

Wavelength (Angstrom) 



Fig. 8. — As Figure [T] but showing the Si xiv Lya line. The vertical dashed lines show the rest wavelengths of the two components of 
the line (6.1804 and 6.1858 A, respectively). 



0.992) the lines grow fainter with respect to the broad- 
band flux. This is what one would expect as the amount 
of absorption starts increasing: the emission lines in the 
~2-3 keV range will be more strongly attenuated than 
the broadband flux over the whole 2-10 keV band. Fur- 
thermore, one would expect the longer wavelength lines 
to show this effect the most. From Figure [12] one can see 
that in CXO030502 this effect is weakest for the S xvi line 
and strongest for the Si XIII resonance line. Rather sur- 
prisingly, however, the Si xiv line is less affected than the 
S XV resonance line. Furthermore, the S xv resonance 
line brightens slightly with respect to the broadband flux 
between CXO030502 and CXOoaoeie- Note in the flnal 
observation, after the recovery (CXO030926; 4> — 1.043), 
that the lines are very faint with respect to the broad- 
band flux. This is because absorption is still having a 
strong effect on the spectrum, and the observed 2-10 keV 
flux is coming from shorter wavelengths than the Si and 
S lines (one can see from Fig. [S] that most of the flux in 
CXOo30926 is shortward of 4 A). 



3.2. Line Shapes 

iBehar et al.l ()2007t l co-added Lya, He-like resonance 
and He-like forbidden lines of Si, S, and Ar and showed 
that the resulting proflle exhibits a signiflcant asymme- 
try on the blueward side of the line. They flnd that 
the lines develop blue wings extending to ~2000 km s^^ 
in CXOo30502 and CXO0306I6, and attribute this to 
the development of a jet with line-of-sight velocity 
~— 2000 km s~^ near periastron. We also looked for ev- 
idence of profile asymmetries, using the individual (i.e., 
non-co-added) lines in each observation. A visual in- 
spection of Figures [7] and [5] suggests that some of the 
lines may indeed be asymmetric. We flnd that some of 
the lines have negative skewness in wavelength (or veloc- 
ity) space, i.e., an extended tail on the blue side of the 
line. For example, the Si Xiv line is skewed in this way 
in CXO021016 and CXO030502 (except in the HEG -Ll 
spectrum). However, the observed asymmetry is not as 
apparent in the Si xiii triplet, which makes it difficult 
to determine whether this apparent asymmetry is real. 
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TABLE 3 
Emission Line Shifts And Widths 



Ion 


Shift Width (FWHM) 




(km s 


(km s ^) 


CXO001119 


S XVI 


-44 ± 76 


1281 ±311 


S XV 


-125 ±42 


815 ±107 


Si XIV 


-34 ±29 


800 ± 82 


Si XIII 


-81 ±23 


555 ± 59 


CXO021016 


S XVI 


-76 ±38 


894 ±127 


S XV 


-77 ±36 


809 ± 89 


Si XIV 


-155 ±24 


1131 ±78 


Si XIII 


-95 ±18 


623 ± 50 


CXOo30502 


S XVI 


-336 ±38 


1333 ±121 


S XV 


-309 ±54 


1316 ±119 


Si XIV 


-514 ±24 


1637 ±68 


Si XIII 


-370 ±27 


1111 ±72 


CXO0306I6 


S XVI 


-387 ±57 


1435 ± 178 


S XV 


-446 ±48 


1281 ±131 


Si XIV 


-665 ±34 


1376 ±92 


Si xiii 


-676 ±41 


1614 ±99 


CXOo30720 


S XVI 


+159 ±260 


1382 ±615 


S XV 


-77 ±89 


821 ± 238 


Si XIV 


-184 ±92 


587 ±257 


Si XIII 


-262 ±59 


894 ± 162 


UXOo30926 


S XVI 


+108 ±127 


507 ± 444 


S XV 


-101 ±71 


708 ± 202 


Si XIV 


-199 ±155 


1213 ±485 


Si XIII 


-122 ±95 


911 ±208 




Fig. 9. — Observed emission line shifts (a) and widths (b) plotted 
against phase. The phase for each observation is taken from Ta- 
ble [T] Phase 0=1 corresponds to the start of the X-ray minimum 
in 2003 June. The gr ay datapoints are for li nes contaminated by 
the CCE component IIHamaguchi et aLll20071) . 




-800 -600 



-400 -200 
Shift (km s'"') 



200 400 



Fig. 10. — The observed emission line widths plotted against 
the line shifts, showing the correlation between the two. The 
ray datapoints are for l ines contaminated by the CCE component 



gray datapoints are tor i n 
IIHa,maguchi et al.|[2007l') . 




33 
X 



O 
> 



Fig. 11. — Observed emission line fluxes plotted against phase. 
The phase for each observation is taken from Table [T] Phase <p = 1 
corresponds to the start of the X-ray minimum in 2003 June. For 
S XV and Si XIII we plot the resonance line flux. The solid line 
show s the 2-10 keV lightcurve measured with RXTE llCorcoranI 
120051) . The gray datapoints are for lines contaminated by the CCE 
component (Hamaguchi ct al. 2007) . 
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Fig. 12. — Observed emission line fluxes divide d by the contem - 
poraneous 2-10 keV flux measured with RXTE l lCorcoranll2005l') . 
For S XV and Si xill we use the resonance line flux. For all lines, the 
flux ratios have been normalized to the values for phase <p = 0.528 
(CXO001119). The gray datapoi nts are for li nes contaminated by 
the CCE component (IHamaguchi et al.ll2007l) . 
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We noted that Gaussians give good fits to the individual 
observed Hnes. A Gaussian profile would give a bad fit 
to a strongly skewed line. 

In order to quantify the amount of asymmetry in the 
observed line profiles, we calculated the skewness of the 
distribution of photon wavelengths that mak e up a given 
observed line. The skewness S is given by (jPress et al.l 
[1991) 

where N is the number of photons, is the wavelengths 
of the ith photon, and A and a\ are the sample mean 
and standard deviation of the wavelengths. If our null 
hypothesis is that the underlying wavelength distribution 
is Gaussian, the standard deviation of S is approximately 
(jPress et al.lll992[ ). In the HETGS spectra, the 
photons are in bins of width 2.5 mA (HEG) and 5 mA 
(MEG). When estimating S, we assume that all the pho- 
tons in a given bin have a wavelength equal to the bin's 
central wavelength. We do not take into account the 
contribution of the continuum, but for most observations 
this should not affect the results too badly, as the lines 
are much brighter than the continuum. 

We looked for skewness in the Lya lines of S xvi and 
Si XIV, and the resonance and forbidden lines of Si xiii. 
We did not include the resonance and forbidden lines of 
S XV, as the S xv intercombination line is more promi- 
nent (see below), which could affect the results. In par- 
ticular we looked for cases where \S\ > 3y^6/N, although 
it should be noted that this might not be a strong enough 
criterion for deciding if the skewness in the line is signif- 
icant We examine the individual HEG and MEG -1-1 
and —1 orders, and also the co- added first-order HEG 
and MEG spectra (to increase the signal-to- noise ratio). 

We found that the Si xiv Lya and Si xiii resonance 
lines are significantly negatively skewed in CXO030502 
and CXOo30616j but the evidence is less convincing for 
the forbidden line in these observations (it is signifi- 
cantly skewed in the MEG -fl order, but not in the 
other orders). The S xvi Lya line is not significantly 
skewed in these observations (see Fig. I13[ which com- 
pares the HEG —1 profiles of S xvi and Si xiv Lya 
from CXOo30502)- For the other observations, there is 
no strong evidence for line skewing - in a few cases a line 
might exhibit skewing in one spectral order, but not in 
the other three. 

Although some of the lines seem to be skewed, visual 
inspection of Figures [7] and [5] suggests that these asym- 
metries are relatively modest. Detailed modeling of these 
line profile asymmetries reveal finer details of the wind- 
wind collision (see iJ5]), but the Gaussian-fitting results 
should provide sufhciently accurate information on the 
gross structure of the wind-wind collision. 

We note that, when comparing the results of fitting in- 
dividual lines from individual orders, the lines are some- 
times slightly offset, possibly due to a slight inaccuracy 
in the position of the zeroth-order image. Also, when 
the lines are analyzed individually, we find that differ- 

IPress et al.l II1992I) caution that "it is good practice to believe 
in skewnesses only when they are several or many times as large as 
[the standard deviation]." 



o 
O 




I ' ■ ' ' ■ 1 

-3000 -2000 -1000 1000 2000 3000 



Shift (km s"'') 

Fig. 13. — Comparison of the HEG —1 profiles of S XVI and 
Si XIV Lya from CXOo30502- For clarity, the profiles have been 
binned by a factor of two, and the S xvi profile has been shifted 
upward by 30 counts. 

ent ions sometimes yield different shifts and widths (see 
Fig. [HI). This suggests that adding the profiles from dif- 
ferent lines and different s pectral orders in o rder to im- 
prove the signal-to-noise (jBehar et al.l l2007[ ) might not 
yield accurate profiles. 

3.3. The TZ Ratios of the He-Like Triplets 

The ratio of the forbidden (/) and intercombination 
(i) line intensities of a helium- like ion, TZ = f /i, can of- 
ten provide useful information on the conditions in and 
location of the emitting plasma. This is because the 
metastable upper level of the forbidden line can be de- 
populated to the upper level of the intercombination line 
by UV photoexcitation or electron collisions: increasing 
the UV flux or the electron density reduces TZ from its 
low-density, low-UV limit TZq. In the case of a hot star 
possessing a stellar wind, where both electron density 
and UV flux vary as TZ < TZq implies that the line- 

emitting region is close to the stellar photosphere. 

In Table [31 we present the TZ ratios for Si Xiii and S xv 
measured from each of our HETGS spectr a. Also in 
the ta ble we present TZq , calculated by Blum enthal et all 
(|l972f ) at the temperature at which the triplet has its 
maximum emissivity (8.9 MK for Si xiii and 14.1 MK 
for S xv), and the UV transition wavelengths to go from 
the upper level of the forbidden line to the upper lev- 
els of the intercombination lines. In all observations, the 
Si xiii TZ ratio is greater than TZq, implying that the 
forbidden line is enhanced with respect to the intercom- 
bination line. This has been observed for O vii in the 
XMM- Newton RGS spectrum of the supernova remnant 
N132D (jBehar et al. 2001), and for several different ions 
in the Ch andra HETGS spec trum of the WR-f-O binary 
WR 140 (|Pollock et al.l[2005h . 

It is possible that too high a continuum level would 
lead to line fluxes that are systematically too low. The 
weak intercombination line would be most severely af- 
fected, and this would lead to an artificially high TZ ra- 
tio. We have investigated whether or not this is the case 
in our analysis by adjusting the range of wavelengths 
we include when fitting to the Si xiii triplet. The re- 
sults in Table [2] were obtained by fitting to the spec- 
tra between 6.4 and 7.0 A (note that the plots in Fig- 
ure [7| do not show this full wavelength range). When 
we use a narrower range of wavelengths, the forbidden 
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TABLE 4 

TZ Ratios for Helium-like Ions 





(A) 








Measured TZ Ratios 








Ion 


(A) 


CXO001119 


CXO021016 


CXOo30502 CXO0306I6 


CXOo30720 


CXOo30926 




S XV 
Si XIII 


738.32 
865.14 


673.40 
814.69 


2.5 ±0.8 
5.0 ± 1.3 


1.2 ±0.3 
6.2± 1.5 


1.7 ±0.4 1.9 ±0.4 
5.7 ±1.5 6.8 ±2.6 


3.1 ± 2.0 
6.9 ±4.5 


1.0 ±0.4 

7.6 ±8.9 


2.0 

2.5 



Note. — TZ — f /i, where / and i are the forbidden and intercombination Une fluxes, respectively. 
^ Wavelengths to go from the upper level of the forbidden line to the upper levels of the intercombination lines - Ai^2 are 
the transition wavelengths for 2^Si 2^Pi,2, respectively (from CHIANTI: IDere et allll997t lYoung et al.ll2003D . 
^ Theoretical low-density, low-UV-flux limit at temperature of maximum emissivity (see eq. [2]; values from 
IBlumenthal et aLllToTl) . 



and intercombination fluxes tend to be smaller. While 
none of the individual decreases is statistically signifi- 
cant, the fact that there is a systematic shift suggests 
that with the narrower wavelength range the line fluxes 
are systematically underestimated. However, we do not 
see the opposite effect when we increase the wavelength 
range from 6.4-7.0 A. The amounts by which the fluxes 
change are much smaller than when we decreased the 
wavelength range, and there is no systematic shift in one 
direction (i.e., some fluxes increase slightly, and some de- 
crease slightly). Furthermore, from a visual inspection of 
the fits, there is no evidence that a power-law is not a 
good flt to the continuum over the range of wavelengths 
that we use. We have also checked whether we can get a 
good fit to the spectra with lower values of TZ by lowering 
the normalization of the continuum model by hand from 
its best-fit value. This should, in principle, increase the 
flux of the weak intercombination line relative to that of 
the stronger forbidden line. However, we find that lower 
continuum levels still lead to values of TZ greater than 
7^0 5 s-iid if the continuum normalization is too low, the 
fit to the continuum is very poor. From these observa- 
tions, we conclude that our large TZ ratios are probably 
not due to an inaccurate continuum level. 

Because of the rather large errors on TZ, only the TZ 
ratios for CXO021016 and CXO030502 differ by more than 
2(7 from TZq, and no observed TZ ratio differs by more 
than 3tT from TZq. However, if we take as a null hypoth- 
esis that TZ = TZq for all six of our observations, this 
gives = 18.35 for 6 degrees of freedom (x^ probability 
= 0.54%). This implies that TZ is significantly differ- 
ent from TZq (as calculated by IBlumenthal et ahl Il972f) 
for at least some of our spectra. However, the combina- 
tion of uncertainties in the ratios and the atomic models 
prevent us from drawing any strong conclusions. One 
might suppose that the large observed values of TZq are 
due to inner-shell ionization of Li-like Si to He-like Si 
(ls^2s + e~ ^ ls2s + 2e~), producing ions in the upper 
level of the forbidden line. This could arise in an ioniz- 
ing plasma, because inner-shell ionization requires both 
a high electron temperature and an abundance of Li-like 
ions, two conditions which ten d not to hold simultane- 
ously in an equilibrium plasma ()Mewe fc Schri ivcr 1978). 
However, Li-like satellite lines would also become impor- 
tant in an ionizing plasma, so that the empirically ob- 
served ratio would no longer reflect just the ratio of the 
He-like ions. Without making detailed NEI calculations, 
it is not safe to make even qualitative predictions for the 
expected TZ ratios of an ionizing plasma. 



For S XV the values of TZ are generally close to TZq. 
In an equilibrium plasma, TZ can be used to place con- 
straints on the electron density ric and the UV flux, and 
hence place constraints on the location of the X-ray- 
emitting plasma. One can express 7?. as a function of 
Tie and the photoexcitation rate (f> to go from the up- 
per level of the forbidden line to the upper level of the 
intercombination line: 



TZ 



TZo 



l + i(P/<j)c) + {n^lncY 



(2) 



where (j)c and nc are quantities dependent only 
on atomic parameters and the electron temperature 
(jBlumenthal et aT]ll972l ). The ratio TZ tends toward the 
limit T Zvi when (jj <^ (pc and rie <C nc- IBlumenthal et ahl 
(|1972f) givenc = 1.9x1014 cm'^ and (?ic = 9.16 xlO^ s'^ 
for S XV at the temperature of maximum emissivity. If we 
assume Mc = 1 0~^ Mp, yr'^ and v^^ 30 00 km s"^ for 
the companion (jPittard fc Corcoranll200^ . we find that 
nc ^ nc everywhere in the companion's wind; unless the 
shock compression ratio is very large (several hundred or 
more), this will also be true in the wind- wind collision re- 
gion. Thus, electron collisions are not expected to affect 
t he TZ ratio. 

IBlumenthal et all ()1972( ) also tabulate 4>i,/4>c, where 
is the photoexcitation rate on the surface of a 10^-K 
blackbody. We estimate 4>-i,/(t>c for 'H Car's companion by 
scaling the Blumcnt hal et al.l (jl972f ) value for S xv, as- 
suming the companion is a blackbody with T ~ 36,000 K 
(this is in the middle of the range of effective tempera- 
tures given by iVerner et all 120051 : 34,000 K < Tgff < 
38,000 K). We obtain 4>*/(t>c = 1-6, which means that on 
the surface of the companion we have TZ — 7^o/(l + l-6) = 
0.77, with TZ increasing toward TZq as we move away from 
the star. Unfortunately, the fact that TZ is fairly large 
even on the surface of the companion, and the large er- 
rors on TZ in Table 21 make it difficult to place strong 
constraints on the location of the emitting plasma. As 
one moves away from the companion, the photoexcita- 
tion rate decreases as (j>{r)/(j)c = 2W{r)(j)i/(j)c, where 
W{r) = 0.5[1 — y/l — (rc/r)"^] is the geometrical dilu- 
tion factor, and r is the distance from the center of 
the companion, whose radius is Vc- Note that rc is 
not well known, though Ushibashi et al.l ([l999) estimate 
rc ~ 5Oi?0. If we take the result for CXO001119, and 
say that the measurements imply TZ > 1.7 (i.e., the la 
lower limit), this gives r > l.6rc for the location of the 
X-ray-emitting plasma. 

The low S XV 7?. ratio for CXO021016 seems to suggest 
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that the S xv emission originates close to the companion 
in that observation. If we were to take TZ < 1.5, this 
would imply r < 1.2rc. However, closer inspection of the 
spectra shows that the S XV intercombination line is no- 
ticeably brighter in the HEG —1 spectrum than in the 
HEG +1 spectrum. This can be seen in Figure [Til which 
shows the S xv triplet from the CXO021016 HEG spec- 
tra, along with the S xv triplet from the CXO030502 HEG 
spectra for comparison. The TZ ratios for CXO021016 from 
the individual HEG orders are 0.6±0. 2 (-1) and 1.9±1.1 
(-1-1), while the TZ ratio obtained from the HEG -fl fit 
results combined with those from the two MEG orders 
is 1.8 ± 0.5. It therefore seems that the low TZ ratio 
for CXO021016 is mainly due to the bright intercombi- 
nation line in the HEG —1 spectrum. We have exam- 
ined the first-order HEG image, using the CIAO tool 
tg_scalejreg to establish the position of the S xv in- 
tercombination line. We find that there is no detector 
feature (such as a hot pixel) or X-ray source which is 
contaminating the intercombination line in the HEG — 1 
spectrum. We have also compared the forbidden and in- 
tercombination line fluxes measured in the -|-1 and —1 
orders of both gratings for each observation. We have 
done this for Si xiii and S xv. In principle, this would 
be a total of 48 comparisons (6 observations x 2 gratings 
X 2 ions X 2 lines). However, as we cannot fit to indi- 
vidual orders in all cases, in practice we find we can only 
make 30 such comparisons. Among these comparisons, 
only the CXO021016 S xv intercombination line measured 
by the HEG differs by more than 2cr between the +1 and 
— 1 orders (the difference is 2.2cr). With the null hypoth- 
esis that the line flux is the same for both orders, the 
probability of such a large difference is 2.8%. Therefore, 
it is not surprising that, among our set of 30 comparisons 
between the positive and negative flrst-order spectra, we 
find one case where the two values differ by 2.2cr. This 
suggests that the large intercombination line flux in the 
CXO021016 HEG —1 spectrum is a statistical fluke, and 
that there is no convincing evidence that the S xv TZ ra- 
tio for this observation really is signiflcantly lower than 
those in the other observations. 

As a final point, it should be noted that the spectral 
type of the companion is not known, and a 36,000-K 
blackbody may poorly represent its UV flux at the wave- 
lengths relevant to the above analysis. A more detailed 
model of its spectrum is required to place more accurate 
constraints on the location of the X-ray-emitting plasma. 

3.4. The G Ratios of the He-Like Triplets, and the 
Ratios of H-like to He-like Lines 

We also measured the G — {f + i)/r line ratios for the 
helium-like Si xiii and S xv triplets. The G ratio de- 
creases with temperature. It is also sensitive to densities 
for n > 10^^ cm~"^, which is well above the range of densi- 
ties expected in the wind- wind collision in rj Car from hy- 
drodynamical simulations (< 10^^ cm"-^; see for example 
iPittard &: CorcoranI [200^ . The measured G ratios are 
given in Tabled and are plotted in Figure [15] (horizontal 
solid lines). Also plotted in this flgure is the temperature 
dependence of the G ratios based fro m the Astrophysica l 
Plasma Emission Database (APED: ISmith et al.|[200ll ). 
version 1.3.1 (curved solid lines). The measured G ratios 
generally imply temperatures of T < 8 x 10^ K for Si xiii 
and T < 13 X 10*^ K for S xv. 




4.95 5 5,05 5,1 

Wavelength (Angstrom) 

Fig. 14. — First-order HEG spectra of the S xv f-i-r triplet 
from CXO021016 {top two panels) and CXO030502 {bottom two pan- 
els). The histograms show the best-fitting Gaussian line models 
obtained by fitting to each spectral order individually. The ver- 
tical dashed lines show the rest wavelengths of the resonance, in- 
tercombination, and forbidden lines (5.0387, 5.0665, and 5.1015 A, 
respectively) . 
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Fig. 15. — G = {(forbidden -I- intercombination) /resonance} line 
fiux ratios and H-like to He-like line fiux ratios for silicon (left) and 
sulfur (right). The horizontal solid lines show the observed values 
of G (Icr limits), while the curved solid lines show the theoretical 
values as a function of temperature from APED. The dashed lines 
show similar information for the flux ratio of the H-like Lya line 
to the He-like resonance line. (For the latter ratio we use the fiux 
of the brighter component of the Lyo doublet.) 
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TABLE 5 
G Ratios for Helium-like Ions 



Ion CXO001I19 CXO021OI6 CXOo30502 CXO03O6I6 CXO03072O CXOo30926 



S XV 0.68 ±0.10 0.74 ±0.09 0.79 ±0.10 0.87 ±0.10 0.94 ±0.29 0.89 ± 0.25 
Si xiii 0.82 ±0.08 0.73 ±0.06 0.75 ±0.06 0.77 ±0.08 0.71 ±0.15 0.97 ±0.34 



Note. — G — {f i) / t , where /, i, and r are the forbidden, intercombination, and resonance 
line fluxes, respectively. 



Also shown in Figure [15] is the ratio of the flux of the 
brighter component of the H-like Lya line to that of the 
He-like resonance line - the horizontal dashed lines show 
the observed values, while the curved dashed lines show 
the theoretical values, also from APED. This ratio is an 
increasing function of temperature, as the ionization bal- 
ance shifts from He-like ions to H-like ions. The observed 
ratios generally imply temperatures of T ~ 12-15 x 10® K 
for silicon and T ~ 20-25 x 10® K for sulfur. 

With the exception of the silicon measurements from 
CXO030720 (which was obtained during the X-ray mini- 
mum, and has much lower signal-to-noise than the other 
observations), the temperatures implied by the G ratios 
are lower than those implied by the H-like to He-like line 
flux ratios. This could be taken to imply that the gas 
is out of equilibrium, as the electron temperature (given 
by the G ratio) is lower than the ionization temperature 
(given by the H-like to He-like flux ratio), suggesting 
that the plasma is overionized. However, an alternative 
explanation is that the H-like and He-like line emission 
originate from different regions of the wind-wind colli- 
sion, with the He-like emission originating from a region 
with a lower temperature than the H-likc emission. This 
is what one would expect from a plasma with a range of 
temperatures (such as a wind-wind collision, where the 
shocked gas near the stagnation point is hotter than gas 
further out). It is also possible that absorption in the 
cool, unshocked winds of the stars is affecting the H-like 
to He-like line flux ratios. As the He-like line from a 
given element is at a lower energy than the H-like line, 
it will be more strongly absorbed by the cool, unshocked 
stellar winds. This would tend to increase the tempera- 
ture inferred from the observed H-like to He-like line flux 
ratio. 

The main conclusion of this section and the previous 
section is that we have not found unambiguous evidence 
of non-equilibrium conditions from the observed line flux 
ratios. However, when we compare the observed line pro- 
flies with those predicted by a model of the wind-wind 
collision, we find that the emitting region is much smaller 
than expected if the wind-wind collision were in equilib- 
rium, suggesting that the wind-wind collision may be out 
of equilibrium. This modeling is described in and the 
results are discussed in i )6.2l 

4. A SIMPLE GEOMETRICAL MODEL OF THE 
COLLIDING WIND REGION 

It is clear from the preceding section that 77 Car's X- 
ray emission lines show variability around the time of the 
X-ray minimum. We first attempt to understand this 
variability in terms of a simple geometrical characteriza- 
tion of the emission region as a conical surface of constant 
opening angle. This analysis has been applied to features 
in optical emission lines from WR 79 to constrain orbital 







Line of sight 












• 











Fig. 16. — Geometrical model of the wind-wind collision in r) Car. 
The solid circles represent the two stars. The cone (with opening 
half-angle 0) represents the wind-wind interaction region (along 
which X-ray-emitting material is streaming at speed v^^). The 
viewing angle 7 is the angle between the line of sight and the line 
of centers. 

and other parameters of the system (lLiihrg|ll997l). and 
also to X-ray emission l ines from WR 140 (|Pollock et al.l 
[2005h and 7^ Velorum (jHenlev et al.l[2005h . 

4.1. Description of the Model 

We assume that the X-ray emission comes from a con- 
ical emission region with opening half-angle f3, whose 
symmetry axis lies along the line of centers with the apex 
pointing toward the primary star, and along which ma- 
terial streams at speed vo- The viewing angle 7 is the 
angle between the line of centers and the line of sight. 
The geometry is illustrated in Figure 1161 Assuming that 
there is no azimuthal velocity component, the centroid 
shift (w) and velocity ran ge (z^max ~ Vmin) of an emis- 
sion hne are given by f|LuhriTl997l : iPollock et "all 120051 : 
IHenlev et al1[2005l) 

w = — uo cos/3cos7, (3) 
"max - Wmin = 2Au = 2uosin/3sin7. (4) 

The viewing angle 7 can be calculated from the orbital 
solution. We first define 5" as the angle between the line 
of centers at the time being considered and the line of 
centers when the companion star is in front; 5" can be 
calculated from the true anomaly <I> and the longitude of 
periastron lu: 

cos^- = cos(<I> + w - 90°). (5) 

If i is the orbital inclination, then 7 is given by 

cos 7 = cos ^I* sin i . (6) 

When comparing the predictions of this model with 
the observed data, v may simply be equated to the shifts 
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TABLE 6 
ri Carinae Orbital Parameters 



Parameter 


Value 


Reference 


To (periastron) (JD) 


2 450 799.792 


(1) 


Period P (d) 


2024 


(1) 


Eccentricity e 


0.90 


(2) 


Longitude of periastron lj (deg) 


275 


(2) 


Inclination i (deg) 


50 


(2) 



Note. — (l) ICorcoraiil II2005D : (2) ICorcoran et al.l l l2001al) . 
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Fig. 17. — The orbit of ?j Car's companion plotted in the rest 
frame of the primary (marked with an asterisk at the origin). The 
positions of the companion at the times of the Chandra HETGS 
observations are marked with black circles. The arrow shows the 
line of sight (L-of-S) projected into the orbital plane. This orbit 
was plotted using the parameters in Table |6] The length scale of 
the orbit is set by assuming masses of 8 0M(t, and SOMg for th e 
primary and the companion, respectively llCorcoran et al.ll2001al) . 

in Table [3) The relation between Av and the measured 
Gaussian line widths is less straightforward. We assume 
that the observed velocities range from ~ (Shift) — 
(FWHM) to w,„ax ~ (Shift) + (FWHM), and proceed 
by simply equating At; in equation Q to the observed 
FWHM. 

The orbital parameters we assume ini tially are given i n 
Table El which are largely based upon [Corcoran et al.l 's 
(l200lj) ana lysis of the RXTE light curve, with a revised 
period from (CorcoranI (|2005i). Note that the time of pe- 
riastron passage Tq in Table [5] is actually the time of the 
start of the X-ray minimum (Corcoran 2005) , which was 
used to calculate the phases in Table [TJ However, as pe- 
riastron is expected to occur near the time of the X-ray 
minimum, assuming the two times are equal has little 
effect on the results. If the time of periastron passage is 
allowed to differ from the time of the start of the X-ray 
minimum, this will only result in the curves calculated 
below being shifted slightly to the left or right. The orbit 
specified by the parameters in Table [6] is shown in Fig- 
ure [171 The length scale of the orbit is set by assuming 
masses of SOM© and 30Mp, for the primary a nd the com- 
panion, respectively ([Corcoran et al.ll2001al) . However, 
the scale of the orbit is not important for our analysis 
- all that matters is how the viewing angle varies with 
time. 

In addition to the orbital elements, we also need to 
assume an opening angle /3 for the wind-wind interac- 
tion region, and a speed vq for material streaming along 
the cone. From hydrodynamical simulations of the wind- 



wind collision in ri Car, we adopt a shock opening half- 
angle /3 = 58° ([Henlevll2005[ ). This is consistent with the 
shock opening angle estimated from the equivalent width 
of the Fe fluorescence li ne measured with XMM-Newton 
([Hamaguchi et aO2007f l. At large distances from the line 
of centers, the velocity along the shock cone vq tends to- 
ward the terminal velocity of the companion star's wind 
(3000 km s~^). However, the observed emission lines 
are likely to origina te from nearer to the line of centers 
([Henlev et al.l [20031 ) - the outer regions are not favored 
for X-ray line emission because (a) the gas number den- 
sity n falls off with distance from the line of centers (and 
the line luminosity scales as rt^) and (b) the gas temper- 
ature also falls off, reducing the populations of H- and 
He-like ions whose lines we are discussing here. However, 
very near the line of centers (where vq is much lower), 
the gas is too hot for most of the observed ions to exist 
in significant amounts, and so the line emission falls off 
here too despite the greatly inc reased density . Usin g the 
line profile model described in Hcnlcv et al. (2003), we 
find that most of the line emission should originate where 
Vq « 2000-3000 km s~\ 

The solid red line in Figure [THl shows the results of 
the geometrical model compared to the observed line 
shifts and widths. The observed variation in the line 
widths is in qualitative agreement with the model in 
that the widths increase around 0=1 and decrease 
again afterward, although the model parameters we have 
used predict larger widths than are observed. However, 
the agreement between the observed and model veloc- 
ity shifts is poor using the model parameters adopted 
above: away from the X-ray minimum, the model pre- 
dicts large blueshifts of '^800 km s^^ , whereas we observe 
much smaller blueshifts of '^100 km s~^, while near X- 
ray minimum, the model predicts redshifted lines, in con- 
trast to the increasing blueshifts which we observe. Some 
of this discrepancy may be due to the assumed values of 
the shock parameters and orbital elements. We consider 
the dependence of the model velocities and widths on the 
parameters vq, (3, i, u and e below. 

4.2. Dependence on the Shock Parameters 

The model line centroids and widths depend on the 
conditions assumed for the boundary surface of the ide- 
alized wind-wind interaction, namely the flow speed wo 
and the cone opening angle (3. Since the flow speed wo 
appears as a multiplicative constant in equations ([3]) and 
([4]), varying wq simply varies the amplitude of the vari- 
ation in the predicted shifts and widths. For example, 
lowering vq by a few hundred km s^^ would bring the 
predicted widths into better agreement with the observed 
widths. However, the discrepancy between the predicted 
and observed shifts would still exist. 

From inspection of equations ([3]) and ((4]), one can see 
that varying (3 will also change the amplitude of the vari- 
ation in the predicted shifts and widths. In particular, 
increasing (3 decreases the amplitude of the velocity shift 
and increases the amplitude of the variation in line width. 
However, if everything else is kept the same, the model 
still predicts redshifted lines around the time of the X-ray 
minimum, instead of the observed blueshifted lines. 

4.3. Dependence on the Orbital Elements 
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Fig. 18. — Comparison of the predictions of the simple geometri- 
cal model (Q with the observed line shifts and widths. The sohd 
red curves were generated from the orbital parameters in Table |6] 
(i = 50°, uj = 275°, e = 0.9). The other red curves show the effect 
of varying the eccentricity: e = 0.85 (short dashed) and e = 0.95 
{long dashed). The blue curves show the effect of adjusting u! to 
l85°, for two different eccentricities: e = 0.9 (solid) and e = 0.95 
(long dashed). The g ray datapoints are for li nes contaminated by 
the CCE component IIHamaguchi et aLll20071) . 

Varying the inclination i varies the amphtude of the 
variation in the viewing angle 7. In an edge-on binary 
{i = 90°), 7 varies from 0° at one conjunction^^, to 90° at 
quadrature, to 180° at the other conjunction, and back 
again to 0°. In contrast, a face-on binary {i — 0°) is 
always observed at 7 = 90°. In general, 7 varies between 
90° — i and 90° -|- i during the course of the orbit. The 
result of this is that varying the inclination i also varies 
the amplitude of the shift and width variations. Max- 
imum variability occurs when i = 90°, and there is no 
variability for i ^ 0° . However, whereas reducing vo or (3 
reduces the predicted widths as well as the amplitude of 
the variation, as i tends to 0° the width tends to vq sin f3 
rather than to zero (see eq. [4]). We find that simply 
varying the inclination cannot bring the model into good 
agreement with the observations. 

We also considered the effect of changing the orbital 
eccentricity. Increasing the eccentricity means that the 
viewing angle changes more rapidly during periastron 
passage. This in turn means that the predicted shifts 
and widths will change more rapidly. As a result, the 
peak at (/) = 1 in the solid red curve in Figure fTSF a) and 
the double-peaked feature at (/> = 1 in the solid red curve 
in Figure fTST b) both become narrower with increasing ec- 
centricity, and broader with decreasing eccentricity. This 
is shown by the dashed red curves in Figure [TH 

Finally, varying the longitude of periastron u has the 
largest effect on determining the phase dependence of 
the velocities in the model. The solid blue curves in Fig- 
ure (TS] show a model with uj — 185°, which means that 
the orbit has been rotated 90° clockwise. In this orienta- 

Depending on which star has the more powerful wind. 



tion the semimajor axis is approximately perpendicular 
to the line of sight, and the companion passes in front 
of the primary just before periastron. One can see that 
this w does yield lines with small shifts away from = 1, 
and with increasing blueshifts as (j) approaches 1. How- 
ever, the increase in the model blueshift occurs too soon 
in phase compared with the observed centroid shifts. In- 
creasing the eccentricity helps by delaying the blueshift 
in phase, and by making the change in centroid veloc- 
ity more rapid near periastron passage. The blue dashed 
curve in Figure [TST a) shows a model in which ui — 185° 
and e = 0.95 instead of 0.9. Although the agreement is 
not formally acceptable, this model is in rough qualita- 
tive agreement with the variation in the line shifts prior 
to the X-ray minimum, though it fails to describe the 
observed variations in line widths. Further adjustment 
of e, uj, and Vq or f3 might further improve the agreement. 

After the X-ray minimum, the new model predicts 
lines redshifted by a few hundred km s^^, whereas 
the observed lines generally have small ('--^100 km s""'^) 
blueshifts. As noted in [J2l the silicon and sulfur lines 
are significantly contaminated by emission from the CCE 
component in the last two Chandra spectra (except for 
S XVI in CXOo30926)- This means that these lines do not 
accurately reflect the centroids of the lines produced by 
the wind-wind collision. However, with this new value 
of Lu the agreement between the predicted and observed 
widths shown in Figure I18r b) is poorer than it was for 
the original model: away from the X-ray minimum the 
new model predicts larger widths than are observed, and 
the predicted widths decrease around = 1, instead of 
increasing. 

In summary, we have shown how adjusting the various 
parameters in our geometrical model for the line shifts 
and widths affects the model predictions. We find that 
by adjusting certain parameters it is possible to bring the 
model into rough qualitative agreement with the obser- 
vations for a subset of the shifts or widths, but we have 
not found a set of parameters which describes both the 
line shifts and variations in line widths simultaneously in 
all of the observations well, though admittedly we have 
not carried out a complete exploration of the whole pa- 
rameter space. However, by seeing how the individual 
model parameters affect the model curves it is not easy 
to see which combination of parameters would bring this 
simple geometrical model into good agreement with the 
observations. 

5. SYNTHETIC LINE PROFILE MODELING 

In the previous section we showed that there is poor 
agreement between the shifts and widths predicted by the 
simple geometrical model, and those that are observed in 
the HETGS spectra of f] Car. With a longitude of peri- 
astron cu w 270°, we can get reasonable agreement with 
the observed variation of the widths, and with u « 180° 
we can get reasonable agreement with the observed vari- 
ation of the shifts. However, we cannot match the vari- 
ation of both simultaneously. Furthermore, when the 
axis of shock cone is nearly perpendicular to the line 
of sight (i.e., 7 w 90°), the above-described model pre- 
dicts broad double-peaked line profiles (with the peaks 
at ~ zLvQSmf]). With the orbital parameters discussed 
above, we expect at least one of our observations to have 
7 w 90°. However, we do not see double-peaked pro- 
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Fig. 19. — The shape of the wind-wind coUision surface in rj Car, 
calculated using the equations of Cant o et al. C1993). The two stars 
are shown by the black circles on the z axis: the primary is at the 
origin, and the companion is at (1,0). Distances are normalized to 
the stellar separation Ducp- The curves to the right show l{x) (see 
eq. [8]) plotted for two different values of Xp^ak {red: Xp^ak = 0.3; 
blue: Xpcak = 0.6). Note that x is actually the distance measured 
along the wind-wind collision surface. The points at x = 0.3 and 
X = 0.6 are marked on the surface with the red and blue circles. 

files in any of our spectra. To address these issues, we 
have developed a more sophisticated model for calculat- 
ing emission line profiles, taking into account both the 
shape of the wind-wind collision region and the varia- 
tion in the speed at whi ch material flows away from the 
stagnation point. iFalceta-Goncalv es ct al. (2006) showed 
that a similar detailed line profile model, including tur- 
bulent broadening and intrinsic absorption was needed 
to fit the phase-dependent, asymmetric C III 5696 A line 
from the WR-l-0 colliding wind binary Br22. 

5.1. Description of the Model 

We calculate the shap e of the wind -wind collision re- 
gion using the results of ICanto et al.l (|T996), who have 
derived equations for the surface of momentum balance 
between two colliding spherical winds. This model is 
for two totally radiative winds with complete mixing be- 
tween them. While this is not expected to be the case 
in ry Car, it provides a useful starting point for model- 
ing the line emission, in particular for determining the 
shape of the surface of momentum balance. We assume 
that the X-ray-emitting region is optically and spatially 
thin, and coincident with the surface of momentum bal- 
ance. The shape of the wind-wind collision surface de- 
pends on the wind momentum ratio i] = M^Vc/ M^iVri, 
and the flow speed along the surface also depends on 
the wind speeds of the stars Vc and Vr^ . For our canon- 
ical model we adopt rj = 0.2, Vc = 3000 km s~^ and 
= 500 km s"^ (jPittard fc CorcoranI [2003 ) . The re- 
sulting shape of the wind-wind collision surface is shown 
in Figure [19 



Using the 

the flow speed along the wind- wind collision surface away 
from the stagnation point tends toward ^--^900 km s~^ 
far from the stagnation point. However, hydrodynam- 
ical simulations suggest that the flow speed in the X- 
ray-emitting region tends toward the wind speed of the 
companion (i.e., ~3000 km s^^). To allow for this, we 



Fig. 20. — Diagram showing the geometry used in the line pro- 
file calculations. Each point along the wind-wind collision surface 
represents a ring of expanding material, comprising material flow- 
ing tangentially along the wind- wind collision surface at speed tit. 
Locally, the flow velocity makes an angle Aocal with the line of 
centers, while the line of sight makes an angle 7 with the line of 
sight. Note that the curvature of the surface has been exaggerated 
for clarity. 

introduce a ve locity s caling factor /„, by which we mul- 
tiply the C anto et al.l flow speeds before calculating the 
line proflle. This scaling factor is a free parameter in the 
fitting described in the following section. 

We assume that the wind-wind interaction surface is 
cylindrically symmetric about the line of centers. There- 
fore, at each point along the wind-wind interaction the 
emission profile is that of an expanding ring. This ring 
of material flows along the wind-wind collision surface 
at speed Vt = fyVcm, where vga e, is the local flow speed 
given by the ICant o et al.l ()1996[ ) equations. Locally, the 
flow velocity makes an angle /Jiocai with the line of cen- 
ters, as illustrated in Figure [20l Note that this is the 
local shock cone opening angle, as opposed to asymp- 
totic value which we used in 21 We assume that each 
infinitesimal portion of this ring emits a Dirac S function 
line profile, shifted according to the line-of-sight velocity 
V. The emission profile e{v) of the whole ring is then 

e{v) oc [v^ sin^ ^locai sin^ 7 - (w + ft cos Aocai cos 7)^] 

where 7 is the viewing angle, defined as before as the 
angle between the line of sight and the line of centers 
(see Figs. [16] and [20]). Note that e{v) goes to infin- 
ity at Wniin = Wt(-sin^iocaisin7 - cos/3iocaicos7) and 
"max = ■!;t(sin/3iocaisin7 - cos/3iocai cos 7) ; e{v) is unde- 
fined outside those velocities. The function e{v) goes to 
infinity because we assume that the intrinsic line profile 
produced at each point on the ring is a Dirac S function. 
In reality, the intrinsic line profile produced at each point 
on the ring will be broadened; we take this into account 
in our calculations by convolving the line profile calcu- 
lated using equation ([7]) with a Gaussian (see below). 



Canto et al.| (jl996D equations, we find that Note also that the integral of e{v) from v- 



to Vn 



is fi- 



nite, and is equal to the line luminosity of the expanding 
ring. 

Using this model we cannot calculate the line emissiv- 
ity at different points along the wind-wind collision sur- 
face self-consistcntly (unlike, say, the X-ray line model 
based upon hdyrodynamical simulations described in 
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iHenlev et al.l 120031) . Instead, we adopt a simple formula 
for calculating the line luminosity as a function of the 
distance x measured along the wind-wind collision sur- 
face from the stagnation point. The line luminosity l{x) 
per unit distance x is given by 



l{x) = 



-ilinc2:2e"("/"- 



(8) 



peak 



where ipcak is the value of x at which l{x) peaks and Liino 
is the total line luminosity, although in this model we are 
only interested in the line shapes, so Lunc is irrelevant. 
The form of equation ([5]) was chosen after some exper- 
imenting wlth_^fiitm£ simple functions to the curves in 
Fig. 2 of IHenlev efall (|2003[ ). The function l{x) encom- 
passes variations in the temperature, the density, and the 
emitting volume per unit x. Some examples of l{x) are 
plotted in Figure [121 Note that mixing with cooler ma- 
terial and/or non-equilibrium ionization may affect the 
form of equation but such effects are beyond the 
scope of the present modeling. 

Our model line profiles are calculated by summing the 
individual expanding-ring profiles from each point along 
the wind-wind collision surface, weighted by the function 
l{x). We convolve this summed profile with a Gaussian 
with FWHM = 100 km s~^ to model thermal broad- 
ening. The resulting profile is then folded through the 
HETGS response for comparison to the observed profiles, 
as described below. 

5.2. Comparison to the Observed Profiles 

The comparison to the observed profiles was carried 
out using XSPEG^^ vll.3.2. We generated a grid of pro- 
files with 7 = 5°, 10°, 15°, 175°, /„ = 1, 1.25, 1.5, 
5, and Xpcak = 0.1, 0.2, 0.4, 6.4. We converted 
the profiles from velocity space to energy space using the 
rest energy of the line we wished to analyze, and used 
the grid of resulting profiles to generate an XSPEG table 
model^^. 

In our analysis we concentrated first on the Si xiv Lya 
line, as the velocity resolution is higher at its wavelength 
than at that of the S xvi Lya line, and there are no prob- 
lems with confusion with nearby lines, unlike the He- like 
f-i-r triplets. Our initial approach was to fit the model 
profiles to the observed lines with 7, /„ and Xpeak all 
as free parameters. We also added a power-law compo- 
nent to model the continuum, and for a given observa- 
tion we fit the model to all four unbinned spectral orders 
(HEG ±1, MEG ±1) simultaneously, using the C statis- 
tic (a modified form of the lCashI [1979 j statistic, which is 
implemented in XSPEG). We applied the model to the 
first four HETGS observations (GXO001119, CXO021016, 
CXOo30502, and CXOosoeie), as the last two (GXO030720 
and GXOo30926) are contaminated by the GGE compo- 
nent as discussed above. 

The best-fitting viewing angles we obtained were sim- 
ilar for all four observations we analyzed: «34° for 
GXO001119, «22° for GXO021016 and GXO030502, and 
R:!l4° for GXOo306i6- This is surprising, given the 
large range of phases over which the observations were 
taken (for example, the phase changed by wO.l between 

http:/ /heasarc. gsfc.nasa.gov/docs/xanadu/xspec/ 
http:/ /xspec.gsfc. nasa.gov/docs/xanadu/xspec/xspccll/manua' 



CXO021016 and CXOo30502, yet the best-fitting viewing 
angles for these two observations differ by Ri0.1°). We 
could not find an orbital solution (specified by to, i, and 
e) which matched the best-fitting viewing angles for all 
four observations. 

We therefore tried a slightly different approach, by try- 
ing to find an orbital solution which would give model 
line profiles consistent with the observed profiles for all 
observations. We fixed e — 0.95, and for a few sample 
values of oj and i we generated theoretical line profiles 
and compared them to the observations, allowing Xpeak 
and /„ to vary until the C statistic was minimized. We 
constrained /„ to be the same for all four observations 
we investigated. Figure W\\ shows these best-fit line pro- 
files for uj = 270° (green) and w = 180° (blue), with 
i = 50° in each case. These values of u and i are si milar 
to the value s publ ished by ICorcoran et al.l (|2001aD and 
ISmith et al.l (|2004D . respectively. They are also similar 
to the values discussed in As shown in Figure [211 
these values of uj and i result in profiles which have too 
much emission redward of the Si xiv line center. This is 
especially true for models in which lo = 180° and i = 50°. 
The ui = 270°, i = 50° models do a reasonable job in 
matching the Si xiv, except for the last observation just 
before the start of the X-ray minimum (CXOosoeie)- We 
then attempted to see if we could generate a reasonable 
fit to all the observed profiles for some value of lu and 
i. After some experimentation, we found that a model 
with UJ = 210° and i = 70° yielded profiles that pro- 
vided reasonable descriptions of the shapes of the Si xiv 
lines in all the observations. These profiles arc shown in 
red in Figure [2T] The orbit of ry Car with e — 0.95 and 
LO = 210° is shown in Figure [22|(cf. Fig. [H]). The best- 
fitting values of and Xp^ak are shown in the upper part 
of Table [7] The Xpoak values imply that the Si xiv emis- 
sion originates further from the stagnation point (relative 
to the stellar separation) in the later two observations, as 
a^pcak is ~8 times larger for these observations. The fact 
that Xpcak is larger just before periastron than at apas- 
tron means that at periastron the Si xiv emission origi- 
nates from a region with much higher fiow speeds than at 
apastron (compare the values of uc96 in Table [7]). This 
explains why the model gives relatively narrow lines for 
CXO001119, even though there is material flowing almost 
directly toward and away from the observer, and why the 
model gives lines blueshifted by a few hundred km 
for GXO0306I6, even though the angle between the flow 
velocity and the line of sight is large (see Fig. [22). 

We repeated the above fitting with the S xvi Lya line. 
In general, it could not discriminate between different 
sets of orbital parameters as strongly as the Si xiv line, 
but of those that we investigated, i — 70°, uj — 210° 
matched the observed S xvi profiles the best. Table [7| 
also shows the best-fitting values of and a;pcak for 
S XVI. The values of for Si xiv and S xvi are in 
good agreement. This is as expected, as fy is a parame- 
ter describing the global flow properties of the wind-wind 
collision, and so it should not be line dependent. 

The best-fit values of imply that material is flow- 
ing along the wind-wind c ollision surface at higher 
speed than is given by the iGanto et al.l (|1996( ) equa- 
tions. Far from the stagnation point, vcge approaches 



-900 km s" 

l/node61.html 



However, the best-fit values of fy im- 
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Fig. 21. — Comparison of our model profiles with the observed Si XIV Ly« lines, from CXO001119, CXO0210161 CXO030502, and CXO0306I6 
(plotted from top to bottom). Each spectral order is plotted separately. The red curves show the profiles calculated with = 210°, i = 70°, 
the green curves the profiles calculated with uj = 270° , i = 50° , and the blue curves the profiles calculated with lo = 180° , i = 50° . 



TABLE 7 

Best-fitting Model Parameters from Synthetic Line Profile Modeling 



i'scp^ VC96^ Si XIV S XVI 

Obs. (AU) (km s-^) 'x~^ 7^^^ ai^pcak'' 1^ 



No absorption - results for uj = 210°, i = 70° 



CXO001119 


29.2 


200 


079+"'""** 
u.ury_o 007 


(2.3 


AU) 


2 29+"'"* 

^■^^-0.02 


U.J^J^y_0.055 


(3.5 


AU) 


acr + U.Ob 
Z.ZO_o j4 


CXO021016 


17.4 


340 


U.14U_o 007 


(2.4 


AU) 




U.iZD_o 013 


(2.2 


AU) 




CXOo30502 


7.00 


820 


0.5741H1? 


(4.0 


AU) 




u.ouz_o 029 


(2.1 


AU) 




CXO03O6I6 


2.73 


850 


0.683tni^ 


(1.9 


AU) 




U.ciDi_o 041 


(1.0 


AU) 





Absorption included in model - results for uj = 270° , i = 50° 



CXO001119 


29.2 


200 


U.O«^_0.008 


(2.4 AU) 


0.01lili;^i5 3.08in? 


n 197+0.U23 
"■^^^'-0.020 


(3.7 AU) 


< 0.021 2.83+0 12 


CXO021016 


17.4 


340 


U.1^4_o 006 


(2.2 AU) 


0.043tHi^ 


n 11 =+0.013 

U.110_o oil 


(2.0 AU) 


no7+o-oi'' 
u.uur_o 007 


CXOo30502 


7.00 


820 


n 14Q+0.007 

U.14ci_o 006 


(1.0 AU) 


102+0 014 

u.iyz_o oil 


127+0012 
U.J-^(_0.oio 


(0.9 AU) 


u.luo_Q 040 


CXO03O6I6 


2.73 


850 


,„7+0.006 
'-'■J^^' -0.003 


(0.3 AU) 


U.yDD_o 038 


(, „ 0+0.013 

U.llii_o oil 


(0.3 AU) 


-in+o ie 

'^•■^^-0.07 



^ Stellar separation in AU, using stellar masses from|C orcoran et al.l Il2001ar . orbital period from|C orcorani | |20Q5|) . and eccentricity e — 0.95 

from ^5?2l 

^ Velocity along the contact discontinuity at a: — a^p^ak for Si XIV, according to lCanto et al.] ( |1996D . 

^ Distance from the stagnation point at which the line luminosity peaks, in units of the stellar separation and (in parentheses) in AU. 
^ Absorption parameter; see equation |9t . 

^ Scaling factor by which the velocities calculated from the lCanto et all (]1996D equations are multiplied before calculating model profiles. For 
each line, the same scaling factor is used for all four observations. 
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CXO0306I6MEG-I 



Fig. 22. — The orbit of r] Car with the eccentricity (e = 0.95) 
from our hne profile modeling. All o ther orbital parameters are 
the same as those used in Figure 1171 The arrows show the lines 
of sight (L-of-S) projected into the orbital plane for two different 
values of the longitude of periastron — that obtai ned w ithout in- 
cluding absorption in the model (oj = 210°; see i|5.2|l . and that 
obtained when th e eff ects of absorption are included in the model 
(oj = 270°; see i|5.3l l. Also shown is a comparison of the stel- 
lar orientations and geometry of the contact discontinuity for two 
phases corresponding to observations CXO001119 (</> = 0.528) and 
CXO0306I6 = 0.992). In each case the red interval shows the 
approximate x range where the line luminosity l{x) (see eq. [8]) is 
greater than half its maximum value for the Si XIV hne. We plot 
results for the version of the model which includes the effects of 
absorption ( i|5.3l lower part of Table [Tjl. The extent of the S XVI 
emission region is similar to that of the Si XIV emission region. 

ply that the speed along the collision surface approaches 
~2000 km s""'^, which is a significant fraction of the 
terminal velocity of the companion (^3000 km s^^; 
iPittard &: Corcoranll2002 ^. It is expected from hydrody- 
namical simulations of the wind-wind collision that the 
flow speed of the shocked gas approaches the terminal 
velocity of the companion far from the stagnation point. 

It should be noted that in the above we do not take into 
account any line-of-sight velocity due to orbital motion. 
The observed velocity profile is actually a combination 
of the projected flow velocity and the orbital velocity 
of the line-emitting region. Orbital motion would make 
the profiles more redshifted if the companion is moving 
away from the observer before periastron, and vice versa. 
However, the fiow velocity dominates, and we find that 
we cannot get a reliable estimate of the orbital motion of 
the line-emitting region from the data. It should also be 
noted that even if we could measure the orbital motion 
of the line-emitting region, it would not be a direct mea- 
sure of the orbital motion of the companion. If we could 
localize the line-emitting region in space, then we could 
in principle relate its motion to that of the companion, 
but in practice this would be difficult to do. 

5.3. The Effect of Bound-Free Absorption on the 
Observed Line Shapes 

iHenlev et al.l ()2003[ ) showed that bound-free absorp- 
tion could have a profound effect on the shapes of the X- 
ray emission lines observed from colliding wind binaries, 
due to differing absorbing columns through the stars' un- 
shocked winds to different parts of the line-emitting re- 
gion. When a colliding wind binary is viewed at quadra- 
ture, in the absence of any absorption the profiles would 
be broad, double peaked, and symmetric about the rest 
wavelength. However, the absorption of the redshifted 
emission from the far side of the system can result in 
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Fig. 



23. — Model Si xiv profiles with (solid) and without 
(dashed) absorption, compared with the MEG —1 Si xiv line from 
CXOo306i6- As in Figure [2T1 the green profiles were calculated 
with uj = 270°, j = 50°, and the blue profiles with w = 180°, 
i = 50°. 



a profile which is strongly positively skewed, with a 
blueshifted peak and a tail extending to the red. 

We have investigated the effect of absorption on our 
profiles by calculating optical depths through the wind 
of the companion star to different points on the wind- 
wind collision. If we assume the wind is spherically sym- 
metric and non-accelerating, these optical depths can be 
calculated analytically (Jgnaca i2001.) . For portions of 
the emitting regions which are viewed through the com- 
panion's wind, we can parameterize the strength of the 
absorption due to the companion's wind with the quan- 
tity 

kMc 



To = 



(9) 



Sep 



where I?sop is the stellar separation and k is the opacity 
at the wavelength of the line of interest. 

We have repeated the fitting described in the previous 
section, with the addition of tq as a free parameter. With 
this modification, we find that we are able to get a good 
fit (judged by eye) to the Si Xiv and S xvi lines for 
CO = 270°, i = 50° and lu = 180°, i = 50° (i.e., we 
do not have to assume a new orientation, as we did in 
the previous section). This is shown in Figure [531 which 
compares the models with and without absorption for 
u = 270° and 180° with the observed MEG -1 Si xiv line 
from CXO0306I6 (we choose this observation to illustrate 
our point because the original model gave a particularly 
poor fit to this observaton for w = 270° and 180°). 

The lower part of Table [7] shows the best-fitting model 
parameters for the Si xiv and S xvi lines with absorp- 
tion included in the model. This set of results is for the 
orbital orientation uj — 270°, i = 50°. Note again that 
the values of for the two lines are in good agreement 
with each other, as expected. Note also that the values of 
a^poak for the two lines are less variable than in the model 
without absorption, and also that they are in good agree- 
ment with each other. The orbit of rj Car with e = 0.95 
and UJ = 270° is shown in Figure [22] (cf. Fig. [TT]). Fig- 
ure [22| also shows the approximate location of the Si Xiv 
emitting region for CXO001119 and CXOoaoeie- 

The results of this modeling are discussed in H6.2[ 

6. DISCUSSION 
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The X-ray emission lines of j] Car show clear variability, 
becoming blueshifted and broader just before the X-ray 
minimum in mid-2003. This variability cannot be de- 
scribed by a simple geometrical model of the wind-wind 
collision in which the emission originates from a coni- 
cal surface with constant opening angle with a longitude 
of periastron of a; « 270°, which is the value consis- 
tent with modeling of t he X-ray 2-10 keV flux variations 
([Corcoran et al.ir2001a[ ) and from analysis of the absorp- 
tion c omponents in He I P-Cygni features (|Nielsen et al.l 
|2007() . However, we found that a more physical model 
which describes the shape of the contact surface and the 
spatial distribution of the X-ray line emissivity along the 
contact surface, and which takes into account absorption 
in the unshocked wind of the companion, was able to 
match the observed HETGS line profiles at all phases 
with uj = 270° and i = 50°. 

Because of the simplifying assumptions inherent in the 
line profile model, and the uncertainty of the input pa- 
rameters, it is possible that other values of lo and i could 
be made to fit the observed X-ray line profiles. It should 
also be emphasized that we did not attempt to find a 
global best-fitting solution, and so uj — 270°, i = 50° 
is not necessarily the best-fitting set of orbital parame- 
ters. Indeed, we found we could also get a good fit to the 
X-ray line profiles with lo — 180°. However, the impor- 
tant result is that we have shown that a colliding wind 
model can explain the observed line profile variations, 
without having to invoke any additional flow component 
(jBehar et al.|[2007f) . 

6.1. Comparison to Line Profiles of Other Massive 

Stars 

Massive stars are believed to produce X-rays via any 
or all of the following processes: by wind-wind collisions 
in binary stars or multiple systems ("colliding wind" 
emission); by intrinsic shocks embedded in the unsta- 
ble, radiatively driven winds ("intrinsic wind" emission); 
and via the magnetic field confinement of the radiatively 
driven wind ("magnetically confined" emission). 

Stars in which intrinsic wind X-ray emission dominates 
the observed emission generally show strong line emission 
but little emission at wavelengths shortward of 3 A. X-ray 
emission from a few of these systems have been observed 
at high spectral resolution. The emission lines are gen- 
erally broad, with HWHM « 1000 km s ~\ which is typ- 
ically half the terminal wind velocity. ILeutenegger et aTl 
(|2006f ) presented a uniform analysis of the helium-like 
lines in four O stars (C Ori, C Pup, l Ori and S Ori) 
and showed that these stars had rather stronger inter- 
combination lines and lower TZ values than we find in 
the 77 Car spectra; typically the TZ ratio was near 2-3 
for the Si xiii triplet, while the TZ ratio was 1.0 ± 0.4 
for S XV in ( Pup (the only star for which S xv could 
be measured). While t Ori and C Ori are binaries, and 
S Ori is a multiple system, none of these four stars show 
any strong evidence of X-ray emission from wind-wind 
collisions, and X-rays from all these stars are believed to 
be dominated by the emission from intrinsic shocks em- 
bedded within the winds of the stars themselves. The TZ 
ratios from these O stars are about a factor of 2 lower 
than the TZ ratios we measured for rj Car (see Table |4]), 
which implies that the minimum radius of the line- 
forming region for these stars is fairly near the stars. 



1.25 < < 1.67 (|Leutenegger et aLll2006D . where 

is the stellar radius. The exception to this is the S xv 7^ 
ratio from CXO021016J which is similar to the value mea- 
sured from C Pup. An intriguing possibility is that we 
are seeing intrinsic emission from the companion's wind. 
If this is the case, it raises the question of why we see 
this effect in only one observation. 

Colliding wind systems generally show thermal X- 
ray emission, sometimes at wavelengths shortward of 
3 A, and may sh ow strong iron K-shell emission. 
iPollock et al.l (|2005[ ) found that in the long-period, ec- 
centric colliding wind binary WR 140 that the intercom- 
bination lines in all the measured He-like lines in that 
star were very weak, and noted differences between the 
weak intercombination lines in WR 140 and the stronger 
i components observed in O star spectra. The intercom- 
bination lines are similarly weak in the sho rter-period ec- 
centr i c colliding wind sy stem 7^ Velorum (jSkinner et al.l 
[200lt IHenlev et al.l[2005l) . 

Like the colliding wind systems, magnetically confined 
winds can show thermal X-ray emission shortward of 3 A. 
The best-studied examples of this class are 9^ Ori C and 
r Sco. Both stars show hard X-ray emission which is 
modulated by the rotation of the star, giving rise to vari- 
ations on much shorter timescales than either the "intrin- 
sic wind" emitters or the colliding wind systems. X-ray 
emission line profiles from shocked gas in "magnetically 
confined" winds are typically fairly narrow and usually 
show symme t ric pr ofiles. An analysis of 9^ Ori C by 
iGagne et al.l ()2005l ) showed that the observed emission 
lines were relatively narrow (~ few hundred km s~^) and 
that the line centroids are close to zero velocity indepen- 
dent of viewing angle. This contrasts with the observed 
variable line profiles we see in 77 Car. 

6.2. Constraining the Interaction Region 

The X-ray line profiles offer the most sensitive diagnos- 
tic of the flow of the shocked gas produced by the wind- 
wind interaction in ry Car, and one of the few spectral 
diagnostics which can be unambiguously localized. We 
have shown that simple geometrical models in which the 
interaction region is a conical surface do not do a good 
job in describing the changes in both line widths and line 
centroids for any assumed orbital orientation (the longi- 
tude of periastron lu), orbital inclination or eccentricity. 
A more physical model based on the Canto et al. (199(|) 
analytical "thin-shell" geometry of the wind-wind col- 
lision interface provides a reasonable description of the 
line proflle shapes at all of the observed phases for val- 
ues of eccentricity which are consistent with analyses of 
the broad-band X-ray fluxes and the He I P Cygni ab- 
sorptions. If we do not include absorption in the model, 
the value of to we derive from the line profile modeling 
(210°) is significantly less than the value oi uj = 275° 
used in the modeling o f the RXTE X-ray lightcurve by 
[Corcoran et al.l (|2001al), although it is close to the value 
of ix) = 200° derived bv llshibashil ()2001 ) from his analysis 
of the RXTE data. However, when we take into account 
the effect of bound- free absorption in the unshocked wind 
of the companion, we find we can get a good fit to the 
profiles with ui = 270°. We use the model parameters 
obtained for this latter model (from Table [7|) in the fol- 
lowing discussion. Our profile analysis is consistent with 
models in which periastron occurs near superior conjunc- 
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tion, i.e., when the companion star is behind rj Car as 
viewed by the observer on earth. This is in contrast to 
the an alysis of ri Car's miUimcter wavelength flux vari- 
ations () Abraham et all 2005a, b ) and the He II 4686 A 
line r adial velocity curve ( Abrah am fc Falceta-GongalvesI 
|2007( ) , which both concluded that periastron occurs near 
inferior conjunction. 

We first discuss the best- fitting values of tq in Table [71 
We can use equation ^ to derive the mass-loss rate of 
the companio n from tq. For the opacity we us e cross - 
sections from iBalucinska-Church fc McCammor 
with a revised He cross-section from lYan et al 



1992 ). 



1991 



and solar abundances () Anders fc Grevessi Il989f) . Us- 
ing the stellar separations in Table [71 and assuming 
Vc = 3000 km s~^ for the win d speed of the compan- 
ion (jPittard fc Corcoraiil I2002D . the best-fitting values 
of To in Table [71 imply mass-loss rates for the com- 
panion of 0.6 X 10-5 Mq yr-i for CXOooiiig to 5.2 x 
10 -5 Mp, yr-i for CXOn.^nfi iR (cf. 10^^ Mq yr-^ derived 
bv lPittard fc Corcoran! 12002 from their modeling of the 
CXOooiiig spectrum, using hydrodynamical models of 
the wind- wind collsion). Rather than indicating a true 
variation in the mass-loss rate of the companion, the ap- 
parent variation may be a result of the emission region 
being smaller at periastron than at apastron. As a result, 
the lines of sight to the emission region are more likely 
to pass through the wind acceleration zone close to the 
companion, where the densities are higher. These higher 
densities would enhance the absorption, giving a higher 
Tq in the fitting, and hence a higher mass-loss rate. 

As the absorption cross-section at the energy of the 
S XVI line is approximately half that at the energy o f 
the Si XIV line (jBalucihska-Church fc McCammonlll992[ ) , 
one would expect the best-fitting values of tq for S xvi 
to be approximately half the corresponding values for 
Si XIV. Given the uncertainties, our best-fitting values 
of To are consistent with this expectation. 

The best-fitting values of Xpeak for Si xiv and S xvi 
are similar to each other, implying they originate from 
similar regions of the wind- wind collision. As Si xiv and 
S XVI form over a wide range of temperatures with signifi- 
cant overlap, this result is not surprising. Indeed, models 
of X-ray line emission in colliding wind binaries based on 
hydrodynamical simulations show that S XVI and Si xiv 
are expected to form in similar regions of the wind-wind 
collision (see Fig. 2 in Henlev et al. 2003). 

We have compared our measured values of a^poak with 
those expected from hydrodynamical simulations of the 
wind-wind collision, wi th stellar wind parameters based 
on those determined bv fPittard fc Corcoranl (|2002f ). Pre- 
vious simulations of the colliding winds in ij Car show 
that dense clumps from the cold post-shock primary 
wind can mix into the hot post-shock secondary wind 
(see Fig. 2 in lPittard fc Corcoranl 1200^ . However, the 
amount of mixing is dependent on the code and resolu- 
tion used, and may be reduced by neglected processes, 
such as magnetic fields. Since the predicted X-ray emis- 
sion is sensitive to the degree of mixing (J. M. Pittard, 
in preparation), we have chosen in this instance to con- 
struct a model which minimizes the mixing betw een the 
winds. We use the line-profile model of Henle v~et al.l 
(|2003f ) to determine the distance from the stagnation 
point at which the X-ray emission is expected to peak. 
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Fig. 24.— The values of 

•^pcak measured from the observed Si XIV 
and S XVI hue profiles (triangles), plotted alongside the values 
of 3^pcak expected from hydrodynamical simulations (circles). In 
all cases the values of Xpcak have been normalized to the stellar 
separation. 

assuming collisional ionization equilibrium. The values 
of Xpcak expected from these models are plotted in Fig- 
ure [24l Away from periastron {(f> = 1), the model val- 
ues of Xpcak are fairly constant. This is as expected 
if radiative cooling is unimportant, because for adia- 
batic colliding wind shocks the structure of the wind- 
wind collision scales self-similarly with binary separation 
([Stevens et al.lll99^ . As a result, we would expect Xpoak 
(in units of the stellar separation, Dgep) to be constant 
with orbital phase for a given line. As the system ap- 
proaches periastron, the increasing post-shock densities 
make radiative cooling more important, and as a result 
the size of the emitting region shrinks. Figure [24l shows 
that the model values of Xpoak rapidly decrease as the 
system approaches periastron. 

Figure [24l also shows the observed values of Xpoak for 
Si XIV and S xvi from the lower part of Table [71 The ob- 
served values of Xpoak do not match the values expected 
from hydrodynamical simulations. Away from perias- 
tron the observed values are ~7-10 times smaller than 
the model values. Also, the observed value of Xp^ak does 
not rapidly fall just before periastron. Instead, Xpeak is 
fairly constant at ~0.1. Therefore, away from periastron 
the observed lines are forming closer to the stagnation 
point than is expected from the hydrodynamical simula- 
tions. This may be because cooling in the shocked gas 
is more rapid than pure radiative cooling, perhaps due 
to inverse Compton cooling or mixing with cool mate- 
rial. This increased cooling would result in a smaller 
emission region than would be expected if only radiative 
cooling were operating. However, this enhanced cooling 
is unlikely to be more important at apastron, when the 
densities are lower, than at periastron - just before peri- 
astron {(j) — 0.992; CXOoaoeie), the observed and model 
values of Xpc^k are in better agreement. Furthermore, sig- 
nificant mixing would lead to velocities along the wind- 
wi nd collision surface approaching the analytical model 
of ICanto et al.l ([1996f ). in which complete mixing is as- 
sumed, whereas we o bserve velocities ~2-3 times larger 
than the lCanto et al.l values. An alternative explanation 
is that the shocked gas is out of equilibrium. The lower 
post-shock densities away from apastron may mean that 
the ionization temperature lags behind the kinetic tem- 
perature after the gas has been shock heated. As a result. 
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the Si XIV and S xvi emission would tend to originate 
closer to the stagnation point than would be expected if 
the gas were in equilibrium - if the gas were in equilib- 
rium, it would be too hot near the stagnation point for 
Si XIV and S xvi. As the system approaches periastron, 
the post-shock densities increase, and the shocked gas 
starts to equilibrate. By itself, this equilibration would 
tend to increase the size of the emission region in units of 
Dsep, as the gas near the stagnation point becomes too 
hot for Si XIV and S xvi. However, as the system ap- 
proaches periastron, radiative cooling becomes more im- 
portant, which tends to reduce the size of the emission 
region in units of Dsep- Our observed values of Xpeak 
suggest that, as the system approaches periastron, the 
increase in the size of the emission region due to the gas 
equilibrating is canceled by the decrease in the size of 
the emission region due to radiative cooling, resulting in 
a roughly constant observed value of Xpcak- More de- 
tailed hydrodynamical modeling, including modeling of 
the post-shock ionization evolution, is required to see if 
the above-described scenario is likely. Such modeling is 
beyond the scope of this paper. However, we can in prin- 
ciple estimate whether or not non-equilibrium ionization 
is likely to be important by comparing the ionization 
timescale to t he flow timescale, using equation (14) from 
iHenlev et al.l ([2005). We find the ratio of the ionization 
timescale to the flow timescale is 0.94 for CXOooiiig, 0.56 
for CXO021016, and 0.23 for CXOo30502- Unfortunately, 
these results are inconclusive: as the ratios are neither 
much greater than nor much less than unity, it is difficult 
to state with certainty whether or not non-equilibrium 
effects are expected in these observations. As radiative 
cooling is more important for CXO0306I6) the appropri- 
ate ratio to calculate is that of the io nization timescale t o 
the radiative timescale (eq. [15] from lHenlev et ani2005D . 
which is 0.01. As this ratio is much less than unity, this 
implies that the system is in equilibrium just before pe- 
riastron, as we have suggested above. 

The fact that Xpeak (in units of Dgcp) is roughly con- 
stant implies that strong cooling is not affecting the 
Si XIV and S xvi emission as the system approaches pe- 
riastron - if it were, the emitting region would decrease 
in size relative to the stellar separation, as regions fur- 
ther from the stagnation poin t would become t o o coo l 
to emit. In contrast to this, iHamaguchi et all ()2007f ) 
and Paper H present evidence of cooling in the X-ray- 
emitting plasma based on excess emission on the low en- 
ergy side of the Fe XXV triplet near periastron, caused 
by rapid cooling driving the hot gas out of equilibrium as 
the stars approach one another. However, it is possible 
that while the high densities near the stagnation point 
result in strong cooling in that region, where the Fe xxv 
emission originates, cooling is less important further out, 
where the S XVi and Si Xiv emission originates. Even if 
the shocked material is cooling rapidly near the stagna- 
tion point, it is still possible to have material hot enough 
to emit S xvi and Si xiv further out, because material 
is being shock-heated all along the shock surface. 

If the shocked gas is in equilibrium, then the innermost 
region, nearest the apex of the shock cone, is too hot to 
emit much silicon or sulfur line emission when there is 
not strong cooling, mixing from clumps of post-shock 
primary wind, or clumps within the secondary wind. We 
take the inner region of the shock cone where the S xvi 



emissivity is less than 50% of its maximum emissivity 
to be the region where the temperature is high enough 
that sulfur is almost completely ionized, and assume that 
the Fe xxv and Fe XXVI emission originates from some- 
where in this region. This assumption constrains the iron 
K-shell region to be on the surface of the wind-wind in- 
terface within about 0.14 AU just before periastron. If, 
as discussed above, the shocked gas is not in equilibrium 
in the earlier observations, we cannot use this method to 
constrain the size of the iron K-shell region near apas- 
tron. 

7. CONCLUSIONS 

We have presented our analysis of resolved silicon and 
sulfur X-ray emission lines from a series of HETGS ob- 
servations of rj Car at key orbital phases. These lines 
originate in the wind-wind collision zone where the slow, 
dense wind of rj Car interacts with the fast, low-density 
wind of a massive companion star. We have shown that 
line profile variations around the orbit are not consistent 
with simple geometrical models of the line forming re- 
gion. A more physically realistic model which takes into 
account the detailed geometry of the contact discontinu- 
ity and allows for variations in the emissivity distribution 
along the shock boundary can produce both the varia- 
tions in the line centroids and the observed changes in 
line width. This analysis allows us to probe directly both 
the temperature distribution along the shock boundary 
and also the flow of the shocked wind of the companion 
away from the stagnation point at the apex of the shock. 

The Si XIV and S xvi emission appear to originate 
from similar regions, which is as expected given the range 
of temperatures at which they are produced. However, 
away from periastron they originate closer to the stagna- 
tion point than is expected from hydrodynamical simula- 
tions. This may be because the wind-wind collision is out 
of equilibrium, and the line emission is originating from a 
region which would be too hot if the wind-wind collision 
were in equilibrium (although it should be noted that we 
did not find unambiguous evidence of non-equilibrium 
conditions from the line flux ratios). Just before peri- 
astron the size of the Si xiv-emitting region is closer to 
that which is expected from a hydrodynamical simula- 
tion, suggesting that the shocked gas has equilibrated 
at the time of CXOosoeie- We find that the flow speed 
along the wind-wind collision surfa ce is ^3x the flow 
speed given by the analytic model of lCanto et al.l (|1996f) . 
This larger flow speed approaches a large fraction of the 
terminal velocity of the companion star's wind far from 
the stagnation point, which is in fair agreement with de- 
tailed numerical hydrodynamic models of the flow, and 
suggests that the mixing of the cold postshock primary 
wind into the hot postshock primary wind is relatively 
unimportant, perhaps due to the presence of magnetic 
fields. 

We can obtain a good fit to the profiles with an incli- 
nation and longitude of periastron similar to those which 
have previously been assumed (i w 50° and uj « 270°), 
although this is only true if we include the effects of ab- 
sorption by the unshocked wind of the companion. Given 
the simplifying assumptions inherent in the model, the 
uncertainty of the various input model parameters, and 
the fact that we did not search for a global best fit, we 
cannot rule out other possible orbital orientations. Nev- 
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crthcless, an important result of this analysis is that it 
shows that colliding wind models can fit the detailed flow 
dynamics as shown by the variations in X-ray line pro- 
files, without recourse to additional flows in the system. 

These results must be considered preliminary since the 
observed line profile variations need to be confirmed as 
dependent on orbital phase rather than simply secular 
changes in the wind. Chandra HETGS observations are 
scheduled for Fall 2008 and late 2008, and ideally more 
will be carried out around the time of the next X-ray min- 
imum, expected in 2009 January. These observations will 
complement ongoing X-ray monitoring with RXTE and 
ground-based optical and radio observations. Meanwhile, 
the line profile model can be improved by the inclusion 
of more realistic absorption from the wind of r] Car and 
from the wind of the companion star, which might give 
further insight into the mass loss rate of the wind from 
the companion, and by more detailed numerical mod- 
els including turbulence to directly determine the depen- 
dence of the theoretical line profiles on the orientation 



and geometry of the colliding wind region. 
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